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Abstract 



We investigate the one-dimensional strongly correlated electron models which 
have the resonating-valence-bond state as the exact ground state. The cor- 
relation functions are evaluated exactly using the transfer matrix method for 
the geometric representations of the valence-bond states. In this method, we 
only treat matrices with small dimensions. This enables us to give analytical 
results. It is shown that the correlation functions decay exponentially with 
distance. The result suggests that there is a finite excitation gap, and that the 
ground state is insulating. Since the corresponding non-interacting systems 
may be insulating or metallic, we can say that the gap originates from strong 
correlation. The persistent currents of the present models are also investigated 
and found to be exactly vanishing. 
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I. INTRODUCTION 



Strongly interacting electron systems have been always one of the most important subjects 
in condensed matter physics. Although rigorous results and exact solutions are useful, they are 
rare. Recently, Brandt and Giesekus Q introduced model of strongly interacting electrons on d- 
dimensional {d > 2) perovskite-like lattices in which the exact ground-state wave functions were 
obtained for a certain range of the parameters. Mielke |Q] showed that the exact ground state 
can be obtained in similar models on a general class of line graphs. Following the line of Brandt 
and Giesekus, models in which the exact ground state can be obtained were constructed by several 
authors These models are conveniently described by the cell construction of Tasaki Q which 

will be reviewed in Sec. ^ (also see Appendix A and Ref. p). Tasaki proved the uniqueness 
of the ground states in this class of models. Not only the ground state but also the singlet-pair 
correlation function in a model on a tree was obtained [^. Bares and Lee Q performed a detailed 
analysis for one of the models of Strack Q . They proved the uniqueness of the ground state and 
exactly evaluated the equal-time correlation functions by a transfer matrix method. 

It was pointed out that the exact ground states have the resonating-valence-bond (RVB) 
structure It is the so-called hopping-dominated RVB states which is different from the 

tunneling-dominated RVB states. The latter have been studied intensively in connection to the 
high Tc superconductivity |llo| , |ll| , |l^ Jl^ Jl^ , p^ , 1 6 1 . Tasaki and Kohmoto |^ studied the difference of 
the mechanism that causes the resonance in the hopping-dominated RVB states and the tunneling- 
dominated RVB states. 

In this paper we shall exactly evaluate the equal-time correlation functions of one-dimensional 
models (Model A, B, and C) which will be defined in Sec. IV. One of the models of Strack which 
was studied by Bares and Lee is called Model B in this paper. We shall use the transfer matrix 
method for the geometric representations of the valence-bond states [13,16,r7|Jl^, which is different 
from that of Bares and Lee. In one dimension, the extension of the formalism to other models which 
will not be included in this paper is cumbersome but essentially straightforward. It is crucial that 
we do not have to treat large matrices as needed in the method of Bares and Lee. (In one of the 
models by Strack, for example, we only treat 3x3 transfer matrices, while Bares and Lee needed 
those of 16x16.) This enables us to give the completely analytical solution with finite amount of 
efforts. 

It is shown that all the correlation functions decay exponentially with distance. The result 
suggests that the existence of a finite excitation gap. It is expected that the excitation gap is 
originated from the structure of the ground state which is described by a collection of local spin 
singlets. The filling factor of the ground state corresponds to that of a band insulating state or 
metallic one in the non-interacting system. The properties of the ground state and the gap are 
completely different from that in the non-interacting system. The existence of the excitation gap 
is expected to be a general feature of this class of models. In a certain range of parameters, Model 
B includes one of the models by Strack [^J^, where we reproduce the results of Bares and Lee. In 
a limit of parameters, it corresponds to a kind of the Kondo lattice regime in the sense that there 
are one localized electron and one conduction electron per a unit cell, where the ground state is 
described by a collection of a local singlet between them. The persistent currents [|£,2C,21,22,23 
are also calculated and turned out to be vanishing. 

The plan of this paper is as follows: In Sec. we review the cell construction of the models. 
In Sec. m we describe the geometric representation of the correlation functions in arbitrary dimen- 
sions. In Sec. IV, we perform a detailed analysis of the one-dimensional models using the method 
in Sec. III. In Sec. |V|, the absence of the persistent currents is shown. Section VI is a summary. 
The reader who is interested only in the physical results may take a look at Sec. |l| and VI for 
general properties of systems in one dimension, and then read Sec. IV A 1 , IV A 6 , IV C 1 , IV C 4 , 
IV D 1, and IV D 4 as examples. The results for the correlation functions are shown in each end of 
subsubsection in Sec. 
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II. CELL CONSTRUCTION AND THE GROUND STATE 



Let us first introduce the solvable models in arbitrary dimension by following the construction 
in Refs. Q and In the present paper, we only consider the translation invariant latticeJll. The 
lattice is constructed from identical cells C„ with n = 1, 2, • • •, A^, where N is the number of the 
cells. The cell Cn is a finite set of sites, where each site r £ Cn (r = 1, 2, • • •, R, wherei R = |C„|) 
is either a site with infinitely large on-site Coulomb repulsion (a C/ = oo site or a d-site) which 
can have at most one electron, or a C/ = site (or a p-site) which can have at most two electrons 
with opposite spins. The full lattice Kj^ [= U^^iCn) is constructed by starting from the lattice 
Ai = Ci, and adding cells C2, Cs, ■ ■ ■, Cn successively. When we add a new cell Cn+i to the lattice 

(= U"^]^Cj), we identify some of sites in C^+i with sites in A„ in a one-to-one manner. We 
note that a cell is not a unit cell in the sense of crystalliography. We denote sites in the full lattice 
Atv by X (x = 1,2, • • •, IAtvI). A site x may belong to several different cells. The correspondence 
between x (g A^t) and r (g Cn) is given by 

x = /(n,r), (2.1) 

where f{n,r) depends on the model under consideration. See Fig. |l| for the correspondence. 
For a cell C„, we associate a cell Hamiltonian 



Hn= Yl ^n^aK^la, (2-2) 

with 



R 

«n,a = E^r"^Cn., (2.3) 
r=l 

where Ar"^ are nonvanishing complex coefficient^ and are chosen independently in each cellS. (In 
Sec. M, we can impose the twisted boundary condition by making use of this property. In Sees. |T| 
and 1^, we only consider cases where all the cells have the same A^. Thus we have translationally 
invariant systems and we drop the suffix n in a1"^ there.) Here c^^ and cj^. are the annihilation 
and the creation operators, respectively, of an electron at site r with spin a =|, j. They satisfy the 
standard anticommutation relations {cl.^^,Cg .^} = 5r,sSa,T and {cj^^o-, c| = {cr,cr,Cs,r} = 0. The 
projection operator which eliminates a double occupancy on d-sites is 

Vn= n (1-^nT^ra), (2-4) 

where Cn;U=oo is the set of C/ = oo sites in Cn and nr.,o- = cj^c^^ is the number operator. This 
represents the infinitely large on-site Coulomb repulsion. The full Hamiltonian is 



"'^We can construct more generalized models whose lattices are not translation invariant. See Ref. [|| for 
such models. 

^Throughout the present paper, \S\ denotes the number of elements in a set S. 

''when a magnetic field is applied, A^"' are complex. The proof of the uniqueness of the ground state in 
Ref. 1^ holds also in this case. 

^When a site x belongs to more than one cells, aI"^ can be chosen independently in each cell. 



9. 



N 



n=l 



We rewrite this Hamiltonian into the "standard form". From the identities [|1|,|5[ 

for r ^ s 

h;-] - nr,i)Vn for r E Cn-U=oo 

for r E C„;[/=o, 

where Cn-u=o is the set of [/ = sites in C„, the cell Hamiltonian (|2.2|) is 



c T' c"!" = — P c"'" c P 



"=1,1 r,s£C„ 



.cr=Tar7^s(eC„) 



+ E E 2|aw| cU,.+ E E hi") 

a=hl r,seCn I 



where 



(Al"))*Ai") forr/s 



An) 
Ar 

Ar 

E 2|a(") 



for r = s and r, s E Cn-u=oo 
for r = s and r, s E C„.[/=o 



From (|2.5|), we have 



where 



H = -Eq -V Y1 E *x,y4,aS,<x^' 



N 



n=l 



ix,y = E for X = /(n, r), y = /(n, s) 

N 

E 



£^1 



n=l 



and 



V=\{Vn. 



n=l 



(2.5) 



(2.6) 



(2.7) 



(2.8) 
(2.9) 

(2.10) 

(2.11) 
(2.12) 

(2.13) 



It was shown in Refs. [yj5| that for the electron number, 2N ^ the unique ground state of the 
Hamiltonian ( |2.5| ) or ( p.lOD has zero energy and is given by 



N 



G.S. 



V n n« 



t 

ra,(T 



, 



(2.14) 



where |0) is the vacuum state. As we will see in ( |3.3| ), we have a single valence bond in each cell. 
Therefore, the filling of the ground state is 'i/Na, where Na is the number of sites in the unit cell. 
In Sees. and 0, we use the second term in ( |2.10| ) as the Hamiltonian and denote it by Hs- The 
unique ground state of the Hamiltonian Hs is given by (|2.14| ) with the energy Eq in ( |2.12[ ). 
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III. GEOMETRIC REPRESENTATION OF THE CORRELATION FUNCTIONS 



We describe the geometric representation of the norm of the ground state and the correlation 
functions which was formulated by Tasaki ||2^ in arbitrary dimensions. The geometric represen- 
tation of the norm is described in Refs. [^,^^ for the lattice composed of d-sites only. Here, we 
have the lattice with both p- and d-sites. The ground state ( p. 14 ) can be written 



N 



N 



n=l r,seCn 



n=l r<sGCn 



, 



(3.1) 



where 



^ / 4,T4a + 4,T4a for r ^ s 



and 



for r 



^ Ac . 



The operator h\ ^ is the creation operator of the valence bond (i.e. a singlet pair) between sites r 
and s if r / s. It creates a doubly occupied site if r = s. They obey the commutation relations 



,A.u\ = 0' 



(3.2) 



where r, s, t, and u are different sites one another. This operator satisfies the relation 6j ^ = 6j ^ 
and the ground state ( ^.l| ) is a hopping-dominated RVB state according to the terminology of 
It is different from the tunneling-dominated RVB states [p, 10, 11, 12, 15, 14, 15,1^. 



Now we rewrite the ground state ( |3.1| ) in a convenient form for diagrammatic evaluations of 
the norm and the correlation functions. The diagrammatic method was first introduced by Rumer 



2q| and Pauling ||2^. We denote a valence bond by {x,y}, and a doubly occupied site by {x,x} 
which is regarded as a self-closed bond (Fig. |2|). Since a self-closed bond is actually a doubly 
occupied site, it is allowed only at p-sites. Let a valence-bond configuration y be a set of N bonds 
constructed by choosing a single bond from each cell. We show examples in Figs, ^(a) and (b). 
The bonds do not share a d-site since it can have at most one electron. A p-site is shared by at 
most two bonds. In this way, the projection V defined by ( |2.4| ) and ( 2.13| ) is automatically taken 
into account. We denote by V the set of all the possible valence-bond configurations. The ground 
state (|3.l| ) is rearranged, and written as 



G.S. 



E n 

V&V {x,y}eV 



A* ftt 



(3.3) 



A. Norm of the ground state 
From (|3.3D , the norm of the ground state is 



($G.S.| '^G.S.] 



vev V'ev \ 



n K',y' n 

{x',y'}£V' {x,y}eV 



bly 



(3.4) 



where \{V) = Yl{x,y}£V ^x,y Let us consider a graph V U V in the expectation value. We call 
a bond which belongs to V as a "bra-bond" and one which belongs to ^ as a "ket-bond" (Fig. 
1^. We only consider graphs VuV in which numbers of "bra-bonds" and "ket-bonds" are equal 
at every sites. (An example is shown in Fig. ^(c).) Otherwise, the expectation value is vanishing, 
since the numbers of the creation and annihilation operators are different at the site. 



The graph V UV can be decomposed into connected subgraphs, since the operators bx'^y' and 
bl.y commute (see (p.2p). An example is shown in Fig. |3|(d)-(g). We denote the number of the 
subgraphs hy n(y L) V'). This decomposition is written 



n(yuv') 
VUV'= UiUUi 



(3.5) 



i=l 



where Ui C V and C// C V. 

By noting that b^'^y' and bl, y commute with each other for distinct x, y, x', and y' (see (pT^)), 
we find that the expectation value in ( |3.4D can be factorized into parts corresponding to connected 
subgraphs. Thus, we have 




n{VUV') j 

n A ([/OA* ([/.)( 
i=l \ 



n 



n 

{x,y}eUi 



(3.6) 



We note that each Ui, in ( |3.6| ) depends on the whole configurations V , V' , and V U V' . 

It sometimes happen that two "bra-bonds" and two "ket-bonds" are connected to a single p-site 
in a connected subgraph Ui U U^. For our calculations, it is convenient to eliminate such sites. This 
is done by using the identity b\. yb]y^^ = —b\^ yb\ ., (see Fig. ^). Examples of eliminations of such 
sites are shown in Figs. H diagrammatically. The procedure in Fig. ||(a) generates a minus sign. 
We assign the sign to the non-closed bond. We shall always apply this procedure hereafter and it 
should be understood implicitly. After this procedure, the subgraph Ui U U[ may be decomposed 
into several graphs. We denote the number of the graphs by n{Ui U U[). The decomposition is 
unique and is written 



n(C/«UC/;) 

u,\ju'i^ J2 Wj u w;. 



(3.7) 



The arrow indicates that Wj, Wj are not necessarily subsets of U U Ul- We only have three kinds 
of graphs Wj U Wj: self-closed bonds (Fig. ^a)), degenerate loops which consist of a pair of bonds 
(Fig. |6|(b)), and non-degenerate loops which consist of even number of distinct bonds (Fig. |6|(c)). 



(3.8) 



We call the graph Wj U W, loop. From the decomposition (p.T]), we have 



n K',y' n biy 

{x',y'}eU' {x,y}<^U, 







n(i/,uc/;) 



where rrij is the number of the procedures shown in Fig. |5|(a) and 



Wi 







n 

{x',y'}eW' 



^x',y' 



n 







(3.9) 



{x,y}<^w, 

We call it weight. The value is classified to 

{1 if Wj U Wj is a pair of the self-closed bonds (Fig. |(a)) 

2 if Wj U Wj is a degenerate loop (Fig. |(b)) (3.10) 

2(— i)'j/2-i if \Yj u Wj is a non-degenerate loop (Fig. |^(c)), 

where Ij is the number of the bonds in the j-th non-degenerate loop. These weights are derived in 
Appendix B. Now the norm of the ground state (|3.6|) is written 



(^G.5.| ^G.S.) = J2 




n((7,UC/;) 

n x{w;)x*{w,)i-ir^w, 



(3.11) 



B. Spin correlation function 



The spin correlation function is given by 

($G.5.|g|Sg|c^G.5.: 



{'^G.S.\^G.. 



(3.12) 



4 \ |)/^ z-component of the spin operator on site x. From (3.3), 



where 

the numerator is written 

{^G.s.\SlSl\^G.s) =Y.T. KV')\*{V)U 
v&v y'ev \ 



n K',J]sisi n &i 



. 



(3.13) 

{m',i)'}6V' / 

From the commutation relations ( |3.2| ) and [5|, hu,v\ = [•S'l, ^1 1;] = where u, and u are different 
sites one another, we can decompose the graph V U V into connected subgraphs Ui U U- as we 
did in subsection [II A| . We only need to consider the case where sites x and y belong to a single 
subgraph denoted as U^^'^'^ U U'^^'^\ Otherwise, the expectation value in ( 3.13| ) is vanishing. The 
decomposition is written 

n{VUV')-l 

VUV = U^'''y^UU'^'''^^ + UiUUi (3.14) 

i=l 

After the elimination of p-sites with four non-closed bonds, we only need to consider the case where 
sites X and y belong to a sing le loop (Fig. 0) denoted as VF(^'2')uVF'^'='^\ Other wise, the expectation 
value in ( p. 13 ) is vanishing. We have 

VUV 

n{u(--v'>UU'^^'y^)-l n(VUV')-l (n{U,VJU[) 



i=l 



\ 



(3.15) 



Mw 



(3.16) 



We note that each Wj, Wj in the second term of ( 3.15| ) depends on the graphs U^^'^^ {j'^^'y^ ^ and 
jji^^y) u jj'i^'V) ^ and those in the third term depends on the graphs Ui, U^, and Ui U U-. We denote 
the total number of the loops Wj U Wj in the second and the third terms of ( p. 15 ) by Mw- Their 
labelings in ( p. 16 ) were rearranged. The weight for the loop W^^'^^^ U W'^^'^^ is 



w{x,y) 



( n K',n.)s'xS'y n bi 



{z',w'}ew 
(-l)'^(^'f)i^O 

4 



n K',^' n t 



(3.17) 



where d{x, y) is the number of the bonds between x and y along the loop and w is given by ( [3.101 ). 
A derivation is shown in Appendix C. From (p.8|), (p. 91), ( 3.16| ), and (3.17), we obtain 

{^G.S.\S'xS'y\^G.S.) 



( -I n{VUV') I 

E (-1)'^(^'^)t n A ([//) A* ([/.) ( 
^,y'ev [ i=i \ 



n K',y' n 

{x',y'}eU', {x,y}eU, 



n{VUV') 



E (-1)^^-4 n 



yy'ev 



n(C/jUC/j) 



(3.18) 
(3.19) 
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where the weight for the loop W^^'"^^ U W'^^''^^ is included in the product. 

C. Correlation function (c^. ,jc|^ ,j) 
We evaluate the correlation function defined by 



From (|3.3|), the numerator is written 



yev y'ev \ 



n K'yjvi'T n 



(3.20) 



O). (3.21) 



We decompose the graph V 'l^V' into connected subgraphs using the commutation relations (|3.2D 
and [c^^CT, = [c^,o-! ^x,?/] — ['i,a^^x,y\ — ^ where z ^ x and z ^ y and obtain ( 3.14| ). After the 
elimination of p-sites with four non-closed bonds, we only need to consider the case where sites 
X and y belong to a single graph and they satisfy one of the following four conditions: (i) site 
x{y) with one non-closed ket(bra)bond; (ii) site x{y) with one non-closed bra(ket)bond and one 
self-closed ket(bra) bond. (See Figs. ^ where a site with a pentagon represents the operator c^^^ 
or cj^^o-.) Otherwise, the expectation value is vanishing. It is convenient to eliminate type (ii) sites 
using the identity cj. ^-ftj. ^ = —c}y^^h\.^^ which is represented diagrammatically in Fig. The results 



are shown in Figs. [lO|(a) and (b). After this procedure, the graph W'^^'''^'^ U W'^^ ' in (3.14) is a 



line where "bra-bond" and "ket-bond" are placed alternately and x' and y' are always at the end 
of the line (Fig. |lO|(c)). The weight is 



w{x',y') 







n 

•{u'.'u'lewf^'j'') 



n 

{u,v}&w('^'y''> 







(_l)'(^'y)/2^ 



(3.22) 



where l{x',y') is the number of the bonds (see Appendix D). Thus we have 



(^G.S.|c^,a4,'xl^G.S.) 

v,V'ev I \ 



n 

■ {z',w'}eu'''''y^ 



^z' ,w' I ^x,cr4,(7 



n 4. 

{z,w}eu('^'y^ 



n{VUV')-l 

n 

1=1 







n K'y n 



^ I A A* (W^-'^y')) (_l)'(-'y)/2+m -Q y |^^^.^ 







J ( ' 



(3.23) 
(3.24) 



v.y&v 



where m is the number of the procedures shown in Fig. H(a) for the subgraph U^^'^^ U U'^^''^\ 



D. Density correlation function 

We first evaluate the expectation value of the number operator 

{^G.S.\n,^,a\^G.S.) 



{nx,a) 

From (13. 31), the numerator is written 



(^>G.s.|^G.. 



(3.25) 
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{'^G.S.K,.\'^G.S.) = E E KV')X*{V)(o 
vgv v'ev \ 



n K'yj'^x,. n 







(3.26) 



• {u',v'}€V' / {«,t)}e\ 
We only need to consider the graph V UV' which contains the site x. Otherwise, the expectation 



value is vanishing. The similar calculation to that in Sec. [Ill B| leads 



n{VUV')-l 

vuV = u u'^""^ + E UiUUl 

4 = 1 

n(!7(^)ut/'(^))-l 



n{VUV')-l 



H/W U VF'^"^ + E WjUW;+ E 



i=l 



E WjUW'j 



(3.27) 



where f/^^) U C/'^''^ is a subgraph with the site x, and W^^^ U VF'^^^ is a loop with the site x after 
the elimination of p-sites with four non-closed bonds. We distinguish two kinds of loops: (i) site x 
with self-closed bonds; (ii) site x with one non-closed "bra-bond" and one non-closed "ket-bond". 
(See Figs. O, where a site with a circle represents the number operator.) (We note that a subgraph 
[/(^) u [/'(^T^ith four non-closed bonds at the site x is decomposed into type (ii) graph and loops.) 
The weight is 



w{x) = h)\ n 

_ J 1 for (i) 

- \ (_i)/{^)/2-i for (ii) 



n 

{z,w}£W('^1 



(3.28) 



where l{x) is the number of the bonds in the graph. A derivation is given in Appendix E. From 
(U), (1^, ( |337| ), and (^^, we obtain 

{^G.S.\nx,cT\^G.S.) 



n{VUV')-l 

- n 



A t/' A* K 



n K',.' n &i 

{u',v'}(^Ul {u,v}(^Ui 



= E A (wf) X* (-i)-^(x) n A (lyj) a*(t^,)(-i)^ 

vyev [ j=i 
The density correlation function is 



{5ncc5ny) = {{n^ - (n^)) (uy - {uy))) 
= {nxriy) - {nx){ny) 

_ {^G.S.\nccny\^G.S.) {^G.S.\nx\^G.S.){^G.S.\ny\^G.S.) 



(3.29) 
(3.30) 



(^G.5.|^G.S.> 
.{<^G.S.\^G.S.) 



{{^G.S.\^G.S.)y 



(3.31) 



where 



P 



D{x,y;a) = -^[^G.s. 



p p p\ p\ _i_ p p -^t 



(3.32) 



From (|3.3|), we have 



D{x,y;a) = Y: E A(0A*(V^) 



X 



n K',v']lU 

■ {u',v'}gV' / V 



c -\- c c 



) n 

/ {u,v}(^V 



. 



(3.33) 



We only need to consider the graph V [JV' which contains the sites x and y. Otherwise, the 
expectation value is vanishing. The similar calculation to that in Sec. IIIB leads ( |3.16| ) where 
is a loop (or two loops) which contain(s) the sites x and y. We classify the loops 
V^^(^'2/) u as (i) sites x and y each belongs to two distinct self-closed bonds, (ii) site x (y) 

belongs to the loop and site y (x) belongs to the self-closed bonds, (iii) sites x and y each belongs 
to two distinct loops, and (iv) sites x and y belong to a single graph. (See Figs. |l2|.) Otherwise, 
the expectation value is vanishing. The weight is 



w{x,y;a) 

Jo 



n 



cM TT 



1 for (i) 

(^_l)'(^')/2-i for (ii) when the site x belongs to the loop 

(_l)'W/2-i X (_i)'.{!/)/2-i for (iii) 

(_l)«(^,s/)/2-i/2 for(iv). 



(3.34) 
(3.35) 



where l{x) and l{x,y) are the numbers of the bonds. A derivation is given in Appendix F. From 
(U), (1^, and (|33B| ), we obtain 

D{x,y;a) 

J2 |a(c/'("'^))a*(c/("'J')) 



x( 



(IT ^z',w' ) 2 ( ^^'"'^V'^'^y^'^^^^'^ ~^ ^^,'^^y, — c^y, — <^^^,cr j Jl^ ^1: 



n(vuy')-i 

X n 



1=1 



X U']X* Ui) 



{z,u>}ec/(= 




(3.36) 



n n 

{u',v'}&Wl {u,v}€Wi 

E < A A* (h^(^'S^)) (-l)™u;(x,y;a) A (p^j) A*(Tyj) (-l)™^'u;j i . (3.37) 



E. Singlet-pair correlation function 

The singlet-pair correlation function is given by 

(^G.5.l4,yVJ^G.S.) 



i^x,y^u,v) 



i^G.S.l'^G.. 



(3.38) 



in 



From ( |3.3| ), we obtain 

l^G.S.) 

= E E Kv')x*{v)(o 

vev v'ev \ 



y{x',y'}eV' 



E 



n(^VUV'U{x,y}U{u,v} 

n A([/nA*([/.)(0 



j=l 



{x,y}ev 



n K',y' biyb^^, n 

^{x',y'}eU' / {x,y}£U 







E 

v,v'ev 



n(^VUV'U{x,y}U{u,v} 

n 



1=1 



n x{w;)x*{w,)i-iy 



(3.39) 
(3.40) 



where V and satisfy the condition that the graph V UV U {x, y} U {u, v} consists of connected 
subgraphs. 
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IV. EXPLICIT CALCULATIONS OF THE CORRELATION FUNCTIONS 



The results of Sees. |n|and III are valid in any dimensions. At the moment, however, a practical 
use of them is limited to models in one dimension in which the transfer matrix method can be ap- 
plied. In this section, the equal-time correlation functions are evaluated exactly for one-dimensional 
models. We shall show the analytical procedures for obtaining them for Model A defined below 
in subsection IV A, illustrating the method in details. The results are shown for a system size 
and in the thermodynamic limit. For Models B and C, we shall only show the results in the 
thermodynamic limit. 



A. Model A 



1. Hamiltonian 



Let us consider a lattice constructed from cells with three sites. We have two lattices which 
satisfy the uniqueness condition of Ref. One of them, which we call Model A (the other is 
called Model B, see subsection IV Q ) is constructed by a cell with two d-sites and one p-site (Fig. 
|l^(a)). Note that a cell is not a unit cell. A unit cell is composed of a d-site and a p-site. In 
the models constructed by the cell construction, Model A is the simplest one for the following two 
reasons. The structure of the lattice is the simplest. (We can construct lattices from cells with 
two sites. The exact ground state, however, contains two electrons per site and is fully-filled.) The 
calculation of the correlation function is easier than that of Model B. 



The cell Hamiltonian ( |2.2| ) is obtained by choosing ai^} = X]r=i ^r-c-. 



+A 



2 "-2,0- 



4 



in (|2.3| ) and setting A3 = 1 without loss of generality (see Fig. |l3|(a) for intra-cell index). Here c\ 
(c^,^) is the annihilation operator on a d{p)-site. The full Hamiltonian is obtained by identifying 
the site 1 in the (n — l)-th cell with the site 2 in the n-th cell (Fig. 0(b)). The Hamiltonian is 



N 



HS = VJ2{J2[ {-X,X2ci+Jci^ - Xici+JcP^, - X2ciJcP,, + h.c.) 

a=1,i n=l 

+eicijci^ + e^<.^c^,, ] + ei^,c%_,J c%^,^^} V, (4.1) 



where the on-site potentials are ef = — 2A2, = — 2(Af -|- A^) (2 < n < N), eff_^i 
= — 1. A unit cell is labeled by n. The ground state is 



-2Af, and 



G.S. 



N 



n 



a 



(A) 



N 

n (Ai<i+A24i,.+<i 



(4.2) 
(4.3) 



which is a half-filled state. 



2. Band structure in the single- electron problem 



Before studying the ground state (|4.3| ), we investigate the corresponding non-interacting system. 
We consider the system with an even number of unit cells under the periodic boundary condition 



-N+l 



of. A one-particle state can be written 



n a=p,d 



(4.4) 
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where are complex coefficients. The single-electron Schrodinger equation H\^f) 
where E is the energy eigenvalue, corresponding to the Hamiltonian (|4.1|) is 



From the Fourier transformation = X^fc ^^^^fk where 

27r 47r f-1 

the Schrodinger equation in the momentum space is 



E^i 



-2(AiA2 cos k + Xl + Xfjifil - (Ai + X2e-'^^)ipl 



E^l = -{X, + X2e^'')ipi-Vl 



The eigen energies are 



E± 



2X1X2 cos A: + 2Af + A^ + 1 =F V (2A1A2 cos A: + 2A2 + 2A|) - 4 (Af + A|) 



E\^f), 
(4.5) 



(4.6) 



(4.7) 



(4.8) 



where — , + are the band index with — (resp. +) corresponding to the + {resp. — ) sign. The energy 



gap between two bands is A = ^{2X1X2 - 1) (1 + 4Af - 2A1A2 + 4A|). When 2A1A2 -1 = 0, the 
gap closes at A; = vr. (See Fig. |l4[) 

In the ground state (|4.3| ), there are 2N electrons. Since there are 2N sites in the lattice, the 
electron number corresponds to full-filling of the lower band. Therefore, the ground state of the 
non-interactiong system is insulating for 2A1A2 — 1 7^ 0. It is metallic when 2A1A2 — 1 = 0. 



3. Norm of the ground state 



Before calculating the correlation functions, we evaluate the norm of the ground state (4.3), 



since the state is not normalized. For the sake of convenience, we draw the lattice shown in Fig. 



lS(b) as Fig. |l3|(c). The ground state admits the geometric representation (|3.3D , where an example 



of the valence-bond configuration V is shown in Fig. p^(a). The geometric representation of the 



norm is ( p. 61 ), where an example of the graph V UV' is shown in Fig. 15(b). We first evaluate 
the contribution from a graph V U V . It can be decomposed into the subgraphs Ui U Ul {i = 1, 
2, • • •, n(y U V')). No loop extends over more than two cells and there is no p-site with four 
bonds, since the sites which is identified in the cell construction are d-sites. We do not need the 
procedures shown in Figs. 0. Therefore, the graph Ui U U- cannot be decomposed further and 



we find n {Ui U U^) = 1 in ( )3.7l ) and mj = in ( ^.11 ). The cells with the graph are classified to 



four kinds (Figs, ^(a)-(d)). Consider the graph shown in Fig. |l|(a). From ( ^loD , the weight for 
the degenerate loop is 2 and the contribution from A ({7^') A* (C/j) in ( p.ll| ) is A2. Therefore, the 
contribution from the graph is 2A2. For other graphs see Figs. |l^. 

The sum over the graph V uV in (^]^) is equivalent to that over all the combination of above 
four kinds of cells under the restriction that a d-site has at most two valence-bonds. (The restriction 
means, for example, that the identification of the right d-site in Fig. p!6|(a) with the left d-site in 
Fig. |l^(c) is forbidden.) Hereafter we shall always take into account the restriction and it should 
be understood implicitly. To evaluate the sum we use the transfer matrix method. We have to 
distinguish two cases due to the restriction. Let An and i?„ be the quantity defined by the right- 
hand side of ( |3.6D on the lattice A„. For A^, the sum is taken over all the combination of the cells 
shown in Fig. |16| with the restriction that the n-th cell is represented by Fig. [l6|(a) or (b). For Bn, 
the sum is taken as was done for A^ with the restriction that the n-th cell is represented by (c) or 
(d). They are represented diagrammatically 



An = (4.9) 

Bn = . (4.10) 
We note that the norm in the system size is 

(^G.5.|«'G.5.) =An + Bn. (4.11) 

Given An-i and B^-i, we can form An and Bn by attaching them to the n-th cell. Let us consider 
An first. When the n-th cell is represented by Fig. p!6|(a), we can attach An-i and cannot attach 
Bn-i- When the n-th cell is represented by Fig. p!6|(b), we can attach An^i or Bn~i to it. Thus we 
have the recursion relation 

An = 

= 2XlAn-l + An-l + Bn^l. (4.12) 

For Bn the similar calculation leads 

Bn = 

= 2XlAn-i + 2\lBn^i + 2Xl\lAn-i. (4.13) 
They are conveniently written in a matrix form as 

(t)-"(t;)' Si .if)- 

The initial vector is / = (^o, ^o)^ = (1, 0)"^, since any cell in Fig. |l6| is allowed as the first cell. To 
obtain the quantity (|4.1lD , we choose the final vector F = (^at, Bp^y = (1, 1)^. Since the transfer 
matrix is not symmetric, it is convenient to diagonalize it using the right and left eigenvectors. The 
matrix can be diagonalized as 

T„ = RDR-^ = L-^DL, (4.15) 



where 




where are the eigenvalues of T„, ei = {2\\ + 2\\ + l + u)i) /2 and 62 = (2Af + 2A2 + 1 — wi)/2 with 



uji = ^4Af + 4A| + 4Af + 4A| + 1. They satisfy ei > 62 > for Ai / and A2 / 0. We choose 

the left eigenvectors Li = (Lii,Li2)^ and L2 = (-^^21, -^22)^ corresponding to the eigenvalues ei 
and 62, respectively, and the right eigenvectors Ri = (-Rii,i?2i)^ and R2 = {R12, R22)^ ■ Using the 
diagonal matrix C = LR, we obtain 

Tn = RDC^^L, (4.17) 

where 

(4.18) 




From the quantities, the norm of the ground state in the system size is 



^G.S. 



=^^T^---T2Ti/ 



F'^RD^C-^L/ (4.19) 
1 + 2Af + 2A^ + 4Af Ai + wi ^ 1 + 2Af + 2Ai + 4Af A^ - uoi 



^2 , 



2uJi 2uji 

where we used the relation C^^LR = I, Here I is the identity matrix. In the thermodynamic limit 
the eigenvalue ei dominates and we have 

\ _ l + 2Af + 2A^ + 4AfAi + a;i 



^G.S. 
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4- Expectation value of the number operator 



We calculate the occupation on a p-site {n\^)- The geometric representation of the expectation 
value is ( |3.29| ), where an example of the graph V U V' is shown in Fig. |l^(c). We do not need 
the procedures shown in Figs. |5|, because there is no site with four bonds. Therefore, we find 
n {Ui U C/j') = 1 in ( 3.27] ) and nij = in ( 3.3[1| ). In the representation, the operator n^^ modifies the 
weight associated with the graph which contains the i-ih cell. Therefore, we replace the transfer 



matrix Tj by N-''^ which is a matrix associated with the operator n^„. The expectation value can 
be written 



G.S. 



G.S. 



F^T^---T,+iN;^)Ti_i---Ti/ 



(4.21) 



We derive the matrix N^^''. We have three kinds of graphs on the i-th cell (Fig. 17 (a)- (c)). Let A^^'^ 
and be the quantity defined by the right-hand side of ( |3.29| ) on the lattice Aj. The restriction 

for the sum is that the i-th. cell is represented either by Fig. |l^(a) or (b) for a\^^ and by (c) for 

(p) 

. They are represented diagrammatically 



A. 



ip) 

i 

(P) 



(4.22) 
(4.23) 



and the recursion relations are 



A 



ip) 



^2^i-l + + Bi- 



B. 



ip) 



From ( 3.28 ), the loop with operator n^.^ has weight 1. It is written 



(p) 

i 

(P) 



N 



(p) 



Bi-i 



N 



(p) 



l + Ai 1 
A? A? 



(4.24) 
(4.25) 

(4.26) 



From (|]51|), ( |452D , ( p^ ), and ( p6| ), we have 



ijy(p)'jii— 1 J* 



C1-C2 


LI 




N-i / \ 


N 


2 1 


■[(l + 2Af)(l + 2A 


D+a;i]-[(l + 2Af)(l + 2Ai)-a;i] (f^) 


1 



(4.27) 



where 



Ci = (ei - 2A?) (ei - 2Xl) (ei - A? - A^) 

C2 = (ei - 2Ai) (ea - 2A?) (dcs - A^ei - Afea) /es 



C3 = (ei - 2A?) (e2 - 2A^) (eies - A?ei - A^ez) /ei 
C4 = (e 



>e2 



2A? 



62 



2X1 



62 



2A2 



2A2 



in the system size N. In the thermodynamic limit, N, i, N 



00, we obtain 



(4.28) 



Since there are two electrons in a unit cell, from (4.28) we have the occupation on a d-site 

We can also obtain the same quantity from the geometric representation of the expectation value. 
From the graphs shown in Fig. 17(d)-(f), we have the matrix associated with the operator nf^ 



N 



id) 



X2 
-^2 







-^1-^2 -^1 



(4.30) 



which will be used to obtain the density correlation function. 

The results are shown in Figs, ^(a) and (b) for a = p and d, respectively. We consider the 
following cases: (i) Ai, A2 > 1 (Ai = A2); (ii) |Ai|, IA2I < 1 (Ai = A2); (iii) |Ai| < 1, A2 > 1. For 
(i), on-site potentials satisfy the relation <^ ep. There is almost one electron per a site. For (ii), 
on-site potentials satisfy the relation < e^, — » A1A2. The d-sites are almost empty. The 
p-sites are almost doubly occupied. For (iii), the system decouples to a collection of the pairs of p- 
and d-sites. 



5. Two-point correlation functions 

We calculate the density correlation function for the p-site. We first evaluate the first term in 
( 13.311 ). The geometric representation of the expectation value is ( |3.36 ), where an example of the 
graph V U V is shown in Fig. 15(d). We do not need the procedures shown in Figs. |5|, because 



there is no site with four bonds. Therefore, we find n {Ui U U^) = 1 in ( |3.15| ) and rrij = in ( 3.37 ) 



In the representation, the operators modify the weight associated with the graph which contains 

(p) 

cells between the cells i and j. From the derivation of the matrix N- , the expectation value is 

D{x, y- a) = F^Tm ■ ■ ■ T.+iN^ ^T.-^i • • • T.+iNS^^T^.i • • • Til. (4.31) 
From ( 4.26| ) and ( 4.31 ), the first term in the right hand side of ( 3.31 ) is 

D{x,y;a) 



G.S.\^G.S. 



t(p) 



{P)r 



nTn-1 ■ ■ ■ T,-+iNf ^T,_i • • • T,+iN,^^^T,_i • • • Ti 



I 



(4.32) 











N~j+i / \ 




- (S) 




4\ 


' [(1 + 2A?) (1 + 2Ai) + u,] - [(1 + 2A?) (1 + 2Ai) - u.,] (^) 







where 



Ci = 


Ui 


-2X1) 




C2 = 


[ei 


-2X1) 


(e2 


C3 = 


(ei 


-2x{) 


[62 


C74 = 


(ei 


-2X1) 


(ei 


C5 = 


(e2 


-2X1) 


[62 



2X1 



xl + xl 



ei 



2A2) (A?ei + A2e2 - 6162) (eie2 - A^ei - Aie2) /(eie2) 



= (ei - 2A?) (e2 - 2Xl) (a? + - ei) (Afei + A^ea - eies) /ei 
= (ei - 2A?) (ea - 2A^) (a? + A^ - 62) (A?ei + A^ea - eiea) /ei 



C8 



62 - 2Aj 



62 



2A2) (a? + A 



62 



in the system size A'^. In the thermodynamic hmit N — j , i ^ 00 keeping \ j — i\ finite, we have 



D{x,y;a) 



(CD 



G.S.l^G.S. 



2 +^ 



2 /e2 V'^'^^'' AfAi [1 + 2(Af + Ai + AK A!)] 

61, 



(4.33) 



Prom (4.28), the first term in the right hand side cancels out with the last term in ( p. 31 ). We 
obtain 



(nr)(np 



62 \ 4Af Aj [1 + 2(Af + Aj + A| + Xj) 



61 



2, ,2 



For the d-site, the similar calculation leads 



(61 + 62)^^1 



(4.34) 



(4.35) 



The density correlation functions take negative value and decay exponentially with distance. We 



show the results in Figs. 19 and 20. For the parameter region (i) identified at the end of subsection 
IV A4| , the density correlation between p-sites is enhanced and that between d-sites is suppressed. 
For (ii) and (iii), they are suppressed. 

Since we have no non-degenerate loop, the spin correlation functions (S^Sj) are vanishing for 
\j — il > 2. Since we have no self-closed bond at the sites where the adjacent cells are identified, 
{bijbk,i) is vanishing for |A; — 2| > 1. 

We evaluate the correlation function {c^^cF^^). The geometric representation of the expectation 
value is ( |3.23| ) , where an example of the graph V ijV is shown in Fig. |l^(e) . We do not need 
the procedures shown in Figs. |5|, because there is no site with four bonds. Therefore, we find 
n (JJi U U[) = 1 in (|3.15| ) and rrij = in ( 3.24 ). In the representation, the operators modify the 
weight associated with the graph(s) which contains the cells between i and j. Let the transfer 
matrix associated with the operator cf^ ('^o-) G^'^^^ (G^"'^^^). We need a new matrix G„ for 
the n-th cell (i -|- 1 < n < j — 1). From these matrices the expectation value can be written 



^-R,(p)rp 



•••Ti/. 



(4.36) 



We derive the matrices. We first consider the matrix G^. It is reduced to a number, because we 
have only one kind of graph (a line) on the n-th cell {i + l<n<j — 1) (Fig. 21(a)). Let Gn be 
the quantity defined by a sum. The sum is taken over VuV on the lattice A„ such that the graph 



consists of loops shown in Fig. 16 on the k-th cell {1 < k < i — 1) and the line shown in Fig. 21(a) 
on the n-th cell {i < n < j). A line with 2n bonds is on n cells, since a cell has two valence bonds. 
The weight for the line is (—1)". We assign —1 to each n cells. In this way, the weight ( |3.22| ) is 
automatically taken into account. They are represented diagrammatically 



Gn 



-AiA2G„_i. 



(4.37) 



We have 



-AiA 



1^2- 



(4.38) 
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For the matrix G^^\ let Cf''^^ be the quantity defined by a sum. The sum is taken over V UV on 



the lattice Aj such that the graph consists of loops shown in Fig. 16 on the k-th cell (1 < k < i — 
and the graph shown in Fig. ^(b) or (c) on the i-th cell. The recursion relation is 



G. 



ip) 



and is written 



(p) 



R,ip) / Aj^ 



.R,(p) 



( -^1 + A1A2, Ai 



(4.39) 



(4.40) 



For G^'^^\ let A^^ and B^^ be the quantity defined by the right-hand side of (|3.23|) on the lattice 
Aj. The sum is taken over the graph V UV' such that the graph consists of loops shown in Fig. |l^ 
on the k-th cell (l<A:<i — 1) and the graph shown in Fig. ^(a) on the l-th. cell {i <l < j)- The 



restriction for the sum is that the j-th cell is represented by either Fig. 21(d) for A!f'^ and (e) for 



B^^\ The recursion relations are 



A 



ip) 



B 



ip) 



AiA2Gj_i, 



and are written 



A 



3 

(p) 



A2 

2 ' 



(4.41) 
(4.42) 

(4.43) 



Prom (|4T9|) , (|^), (|438|) , (glo|) , and ( psj) , we obtain 



pTrj.N-jQL,{p)Qj 



-lQR,{p)r^i-lj 



pTr^Nj 



{AAA) 



i-XiXit 



Ci 




-C2 


Is) 


j-1 






N 


LVl < 


■[(l + 2Af) (l + 2Ai)+^i] 


- [(1 + 2A2) (1 + 2A2) 


l(ff) 





where 



Ci = 




C2 = 




C3 = 









2A2) (ei - A?) (ei - 2A2) (ei - A^) /ei 



A?) (ei - 2Ai) (ei - A^) /, 



'1 



in the system size A^. In the thermodynamic limit N — j, i ^ 00 keeping | j — i\ finite, we obtain 

Ai A2 \ (2A2 + 1 + wi)(2Ai + I + L01) 



ei 



AeiuJi 



(4.45) 



18 



Next, we evaluate the correlation function (cf^^c^^^). We show an example of the graph V UV' 
in Fig. |l^(f). (We have only one kind of the line.) Using the graph of the ends of the line (Figs. 
21(f) and (g)) the matrices in (|4.36| ) are 

Gf ('^^ = ( -A,A„ O), Gj'^'^=(\X (4.46) 



From these matrices we obtain 



uji V ei 



The correlation functions decay exponentially with the oscillating sign. 

For a finite lattice under open boundary condition, the system is not translational invariant. 
In the thermodynamic limit, however, by the Fourier transformation of the correlation function 
(Cj^^Cj we obtain the momentum distribution function for a = p, d 



where /q"^ = {nf^) and 



K.> = /"^mO + /i"\ (4.48) 



, 2r cos A; - r 

F(k,r) = ' (4.49) 

^ ' ^ l + r2-2rcosA;' ^ ^ 



where /(p) = -(2Af + 1 + a;i)(2Ai + 1 + L0i)/AeiL0i, /('^^ = l/cui, and r = -AiA2/ei. The results 
are shown in Figs. |2^ and ^ There is no singularity in {n%„)- For the parameter range (i) in 
subsection [V A4 , the momentum distribution for the d-site is completely flat, while that for the 
p-site has a broad peak around A: = 0. For (ii), the momentum distribution for the p-site is almost 
unity for every k and is completely flat, while that for the d-site is almost zero. For (iii), both of 
them are completely flat. 



6. Discussion 



All the correlation functions under consideration decay exponentially with distance. These 

results suggest the existence of a finite excitation gap. Therefore, it is expected that the state 

does not exhibit a metallic state but rather an insulating one. The correlation lengths are given by 

1 

for the correlation 



In 



Inf^ 



A2 



for the density correlation functions and 

functions (c"^c^^) (Fig. The correlation lengths of the correlation between p-sites and that 
between d-sites are the same. The spin correlation functions vanish for |j — i\ > 2. The singlet-pair 
correlation functions vanish for any |j — i\. 

We consider the limit |Ai|, IA2I ^ 1 (Ai = A2). We obtain ^cc = 0. The density correlation 
functions for the nearest neighbor sites vanish. The ground state is described by a collection of the 
decoupled p-sites which are doubly occupied. 

We consider the limit Ai, A2 ^ 1 (Ai = A2). The correlation length converges to a finite value: 



ln( 



2+^2 



and 



ln( 



4+2v^ 



-1 



The density correlation function for the nearest 
neighbor p-sites remains finite, while that for the nearest neighbor d-sites vanish. Since there is 
almost one electron per a site, the correlation between d-sites is suppressed and that between 
p-sites is enhanced. 

For |Ai| <C 1, A2 3> 1, the correlations are suppressed. This is due to that the system decouples 
to a collection of the pairs of p- and d-sites. 

The ground state ( |4.3| ) is a half-filling state. In the non-interacting system, the filling factor 
corresponds to that of a metallic state at 2Ai A2— 1 = and that of a band insulator for 2A1A2 — 1 7^ 0. 



1Q 



Therefore, we have a metal-insulator transition when Ai, A2 are fixed to satisfy 2A1A2 — 1 = 0, and 
the on-site Coulomb interaction U on d-sites is varied from to cxd. In the non- interacting system, 
the correlation length, which is proportional to the inverse of the energy gap, takes a finite value 
for |Ai|, IA2I <C 1, it diverges when 2A1A2 — 1 = hold, and it goes to zero for Ai, A2 3> 1. These 
properties are completely different from that of the ground state ([4.3|). 



B. Notation of the transfer matrices 

In order to present the calculations in latter models efficiently we fix some notations. From the 
derivation in subsection ([VA3), the norm of the ground state can be generally written 



i^G.s.l^G.s.) =F^Tn---T2TiI, (4.50) 

where the matrices depend on the model under consideration. 

We describe the expectation value of a local operator Oi by using the transfer matrices. In 
the geometric representation, when there is an operator Oi, the weight associated with the graph 
which contains the i-th cell is modified. Therefore, we replace the transfer matrix Tj by Oi which 
is a matrix associated with the operator Oi. The expectation value is written 

($G.s.|0*|^a5.) =i^^T,v-T,+iO,T,_i---Ti/. (4.51) 

When Oi is the number operator nf^ where a = d[p) for a p((i)-site, let the corresponding matrix 
be 



Oi = nr. (4.52) 



r(") 

i 

For the two-point correlation function (OfOj), the weight associated with the graph which 
contains the cells between i and j is modified. Let be the transfer matrix between the sites i + 1 
and j — I, and Of (Oj) be the matrix associated with the operator Of (Oj")- We have 

{^G.s.\O^Of\^G.s.) = F'^TnTn-i • • • T.+iOj-P,-! • • • Pi+iOf Ti_i • • • tJ. (4.53) 
We use the following notations: 





= Tfc, 


Of 


= N^"^ 


Of 


= 


Pfc 


= Sfc, 


Of 


— ' 


Of 


j 


Pfc 


= Hfc, 


Of 




Of 


Pfc 


= Gfc, 


of 




Of 





for the quantity D{i,j;a) in ( 3.32 ) 

for the spin correlation function 

for the singlet-pair correlation function 



(4.54) 



for the correlation function (Cj^cj^). 
For the singlet-pair correlation function (&Jj&fc/), we distinguish the following four cases by 7: 

(i) 7=p sites i and j are p-sites and i = j; 

(ii) 'y=pp sites i and j are p-sites and i ^ j; (a 

(iii) 'j=pd site i is p-site and site j is d-site; 
^ (iv) '^=dd sites i and j are d-sites. 

We can evaluate multi-point correlation functions for operators which are constructed from 
fermion operators. The numerator of the correlation function is obtained by an insertion of the 
transfer matrices which are associated with the operators. 
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C. Model B 



1. Hamiltonian 

Model B is constructed from a cell with two p-sites and one d-site (Fig. p^(a)). This is one of 



the models of Strack [Q] which was studied by Bares and Lee [Q. The cell Hamiltonian (2.2) is 

obtained by choosing al^J = J2r=i \(^r,c! = o- + -^2 a- + -^3 '^3,0- (2^) ^'^'^ setting A3 = 1 
without loss of generality (see Fig. |2^(a) for intra-cell index) . The full Hamiltonian is obtained by 
identifying site 1 in the (n — l)-th cell with site 2 in the n-th cell (Fig. ^(b)). The Hamiltonian is 



AT 



u=^,l n=i 



+ 4^+1 



where the on-site potentials are = —X^, = — (Af + A2) (2 < n < N), e^^-^ = — Af, and e'^ = —2. 
A unit cell is labeled by n. The ground state is 



N 



B 

G.S. 



vu n « 



n=l (T=T,J, 
N 

n (Ai<i + A24:i,. + <i 



, 



(4.57) 



which is a half-filled state. In the parameter space, A2 = — Ai, this model reduces one of the models 
in Ref. Q. The model in Ref. |Q is recovered by setting Ai = — A2 = V~^. 



2. Band structure in the single- electron problem 



We investigate the single-electron problem for the Hamiltonian (4.56). We consider the system 
with an even number of unit cells under the periodic boundary condition. The similar calculation 
to that in Sec. [V A 2| leads to the dispersion relations 



2A1A2 cos A; + A? + Ai + 2 =F ^J (2A1A2 cos A; + Af + A^) V 4 



(4.58) 



where — , + are the band index with — {resp. +) corresponding to the + (resp. — ) sign, and k 
is the wave vector with (|4.6| ) . The energy gap between two bands is A = ^ ^/ (Ai + A2)^ + 4 + 
^\/(Ai — A2)^ + 4 — 2A1A2, which is nonvanishing for any finite Ai, A2. 

The electron number in the ground state ( 4.57] ) corresponds to full- filling of the lower band. 
Therefore, the ground state of the non- interacting system is insulating. 



3. Correlation functions 
From ( ^.501 ), the norm of the ground state is 



^G.S. 



^G.S. 



— (i?llLii + 2i?2iL2l + R31L31), 
Cl 



(4.59) 



where the corresponding matrices in ( [4.14 ), ( 4.16| ), and ( [4.18| ) are 
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D 



/ 2Ai + A| 
Af + Af A2 




2Ai 

2Af + Af a: 

2A| 



and C 




F 




(4.60) 



The matrix T„ and the initial and the final vectors are derived in Appendix G 1. Here Cj [i 
and 3) are the eigenvalues of T„ 



1, 2, 
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1 




- tp3^ 




ei = 


3 + 


3 


pi - 




62 = 


s 

3 ~ 


1 

6 






+ i\/3 (p3 + tp3^ 


es = 


s 

3 ~ 


1 

6 


(pi 




- i\f?> [ph + tp3^ 



where s 
Af A2 + A1A2 



2Ai + 2A2 + Af + A1A2 + Ai, t = -4Af - 2A1A2 - 4Ai - 4A? - 2Af A2 - 2AiAi - 4Ai - A| + 



A|, and p = (pi + 32^^)72. Here pi = 16Af + 12Af A2 + 12AiAi + 16Ai + 24Af + 
18Af A2 + 6Af Ai + I8A1 Ai + 24A| + 12Af - 12A? A^ - 12A? A^ + 12Ai + 2\\ - 3Af A2 + 14A? Af - 3Ai A| + 2A^ 



-16A? 



16A2 



2-48Af-32AfA2-32AiAi- 



- 48Ai - 68Af - 72Af A2 - 60Af A^ - 72AiAi ■ 



and p2 = Af Aiv-J-u/X]^ 

68A| - 56Af - 40Af A2 - 32Af A^ - 32Af Af - 40AiA| - 56Ai - 28Af + 12Af A2 + 12A|Ai + 8Af A^ + 
12Af A| + 12AiAi - 28A^ - 8A[ + IGAf A2 - 8A|Ai - 8Af A| + 16AiA| - 8A^ - A? + 4AlA2 - 4A^Ai - 
4Af Ai + lOAf A| - 4A?Ai - 4Af A^ + 4AiAi - \\). They satisfy ei > 62 > eg > for Ai / and 
A2 7^ 0. The matrix L = (Lij) (R = (Rij)) is constructed from the left(right) eigenvectors. We 
choose the left eigenvectors Li = (Ln, L12, ^13)"^, L2 = (L21, L22, -^^23)^, and L3 = (L31, L32, -Lss)'^ 



Af(Af-e,)(l + Ai) 

Li2 = (Af — e,)(2A9 — At — e,), and L^g = Af (2Ai — A2 — ej). We choose the right eigenvectors R 



corresponding to the eigenvalues ei, 62, and 63, respectively, where Lji 

^32 



{Ru, R21, Rsi)^ , R2 



(i?i2, i?22, -^32)^, and ^3 



R2j = (Af - ej){ej 



2Ai 



A2), and R, 



3j 



2A|(Ai + A2 



(i?i3,i?23,^33)^, where Rij 



2Ai(At-e,), 



We evaluate the expectation value of the number operator (nf^). The transfer matrices asso- 



ciated with nf^ are 





^Ai + A! 


2Ai 






1 


^2 


Ai 




2 A1A2 


Af + AfAi 


i) 




2-^1 


A? 








Af 






\ 











(4.61) 



They are derived in Appendix^ From ( |4.27D , (|]5|), ( |4.6q ), and (|46l| ), we obtain 

^j^a ^ ^ I (-^11 ~ ^AiLi2) A2-R11 + {L\2 - F\z)R2\ + Li3i?3i for a = p; 

1 (A2-^^ii + ^AiLi2)i?ii + (A2-L11 + AiLi2)i?2i + ^AiLi2i?3i for a = d, 



(4.62) 



in the thermodynamic limit, N , i, N — i ^ 00. It can be verified that there are two electrons per a 
unit cell as (n^) + {nf) = 2. The results are shown in Figs. 27. The case A2 = |Ai| was disscussed 



in Ref. Q. We consider the case |Ai| <C 1, A2 ^ 1. The occupation on p- (d-) site have a minimum 
(maximum) at a intermediate value of Ai. 

From ( |33l| ), ( |427D , (p^ ), ( ^391 ), (|i!60| ), and ( ^^l] ), the density correlation functions are 



eie2CiC2 



r (a) 
-fi-2 ^2 



+ 1^ 

ei 



7->(") r(") 
^3 ' 



(4.63) 



where 
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rip) 



(Lii — L12/2) Rij + (L12 — Li3)i?2j + LisRsj 



RuLji + {-Ru/2 + R21) Lj2 + {-R21 + RzijLjz 

{\2L11 + A1L12/2) Rij + (A2ill + A1L12) i?2j + AiLi2-R3j/2 
A2(-Rll + R2l)L21 + Ai (iill/2 + i?21 + /2-R31) -^^22, 



in the thermodynamic hmit N — j, i ^ 00 keeping \ j — i\ finite. The density correlation functions 



take negative value and decay exponentially with distance. We show the results in Figs. 28 and 



2£. For I All <C 1, A2 ^ 1, the density correlations are suppressed. They have a minimum at a 
intermediate value of Ai. 

We evaluate the spin correlation function. The transfer matrices are 



(lAfAl, 0, 0), 



„R,{d) _ / 1 ^ 2 1,2 1\2\ 



gi,(p) 



( ° 

V- 



At 



(4.64) 





V 



They are derived in Appendix P 3|. From ( |4.53[ ), ( [4.54| ), ( |4.59D , and ( |4.64| ), we have 



pT'YN-i^L,{a) gj-i-lgH,(/3)rpj-lj 



pTr^Nj 



ei 



(L12 + Af L13) 



X < 



2e|Ci 

'2eTc7-^ll + 2i?21 + Ru] 

\2 



for (3 = p and a 
for (3 = d and a 



2ASeici 



(L12 + AfLia) {Rii + 2i?2i + -R31) for /? = p and a = d 

for (3 = d and a = p. 



(4.65) 



(4.66) 



-2^(^11^11) 

We note that {S^'^^Sp'^) / {Sp^Sp'^)^ because the Hamiltonian is not invariant under the reflection 
of the lattice. All the spin correlation functions take negative value and decay exponentially with 
distance. We show the results in Figs. For |Ai| ^ 1, A2 ^ 1, the spin correlations are suppressed. 

We evaluate the singlet-pair correlation function ( 3.3^ ) where i and j {k and /) are in the same 
cell. For ( [4.55D -(ii), the correlation function is vanishing. We evaluate them for ( |4.55D -(i) and (iii). 
(The position of the operators are shown in Fig. ^(g) and (h) in Appendix |G 4| , where a double 
solid (broken) line represents the operator blj (fefc,/)-) The transfer matrices are 

H„, = 2 



H 



R,{pp) 



1, 0, 



H 



L,{pp) 



H 



R,{dp) 



2A?A2 ( 1, 1, 



H 



L,{dp) 



2X1 
2X\ 



2Ai A2 




(4.67) 



They are derived in Appendix |g1 From (|453|) , (|]5|), (|459|) , and (|46^ ), we obtain 



^TrpAr-A;jji,(7)jjfc-i-ljj-R:(7)rj,i-lj 



pTr^Nj 



^1^2 

ei 



-\i-k\ 



LrsRn 



for 7 



4A1A2 

eici 



Li2{Rii + R21) for 7 



P 

dp. 



(4.68) 
(4.69) 



The singlet-pair correlation function decay exponentially with distance. The results are shown in 
Fig. 31 for 7 = p. For |Ai| <C 1, A2 3> 1, the singlet-pair correlations are suppressed. 
We evaluate the correlation function (cf^cjj^). The transfer matrices are 



-A1A2 



,R,{d) 



-Ai 



'n 

1 + 1 



AiA 



1^2 



A? 



A? 



Ai Ai 
Af 2X1 Af 



1 + Ai 1 
A? A? 



,L,{d) 





Ai 


A? 


+ ^AfAi Af 




Af Af 




Mi \ 


A2 


1\2 
2-^1 2^1 




[ ) 



(4.70) 



They are derived in Appendix |G 5| . The matrix G„ is diagonalized 

Gn ~ — A1A2RDC ^L, 



(4.71) 



where 



D 



51 
92 



R 



Rii R12 
R21 R22 



Lii L12 
L21 L22 



ci 
£2 



(4.72) 



Here gk are the eigenvalues of G„, gi = {1 + Xi + X2 + W2)/2 and 52 = (1 + Ai + A2 — uj2)/2 with 
0J2 = y/1 + 2(Ai + A2^+ (Ai - A2^2. They satisfy gi>g2>0 for Ai / and A2 / 0. We choose 
the left eigenvectors Li = {Lii,Li2)^ and L2 = (^^21, ^22)"^ corresponding to the eigenvalues gi 
and (72, respectively, where Lji = gj — Af, Lj2 = 1, and the right eigenvectors Ri = (i?ii,i?2i)^ 
and i?2 = {R12, -^22)^ where = gj — Af , i?2j = Af . Here C = LR. 
From these matrices and ( |4.44 ), we obtain 



-AiA 



1^2- 
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'ei 



-AiA 



1^2- 



Ac) 

J2 ' 



(4.73) 



where 



p{p) 



J m 



CiC„ 



\{9m 



Af )Li2 + Af-Li3 



(5„ - Xl)Rii + R21 



1 



eigmCic„ 



X2{9m — Af)Lii -I- -grnXfLu 



Xligm - Af ) (i^ii + R21) + Af (i?ii + 2i?2i + R31) 



(4.74) 



(4.75) 



The correlation functions decay exponentially. 

By the Fourier transformation of the correlation function (q^c^^), we obtain the momentum 
distribution function for a = p and d 



(n?) = /WFi(fc,n) +/ri^2(fc,r2) +/, 



' 



(4.76) 



where /q"^ = {nf^) and 



2rj[cos k — ri] 
1 + rf — 2ri cos k 



(4.77) 



where rj = —XiX2gi/ei- The results are shown in Figs. ^ and 33. There is no singularity in {n'^ „)■ 
For |Ai| <C 1, A2 3> 1, the momentum distributions for the p- and d-sites are flat. 
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4- Discussion 



All the correlation functions under consideration decay exponentially with distance. These 
resluts suggest the existence of a finite excitation gap. The existence of the energy gap is numerically 
confirmed in Ref. |2^. Therefore, it is expected that the state is insulating. The correlation lengths 

are given by S,nn = ( ff ) , Css = ^bb = In ( ) i ^'^^ the spin and singlet-pair correlation 

functions, respectively, and ^cc = (Fig- H)- (They satisfy the relation ^cc > ^nn 

> S^ss = S,bb-) The correlation lengths of the correlation between p-sites and that between d-sites 
are same. We note that the spin correlation is ferromagnetic. 

We consider the region Ai, A2 ^ 1- The d-sites are almost empty and the p-sites are almost 
doubly occupied. The correlation lengths behave ^cc ~ Ai and S,nn, (,ss, S,bb ~ ^ (for -^1 = ^2). The 
correlation functions for the nearest neighbor sites are suppressed. 

For I All, IA2I ^ 1, there is almost one electron per a site. The density correlation function for 
the nearest neighbor p-sites are enhanced, while that for d-sites are suppressed. The spin correlation 
function for the nearest neighbor p-sites are suppressed, and that for d-sites are enhanced. There- 
fore, the electrons on the d-sites have a tendency to behave like a localized spin. It corresponds 
to a kind of Kondo lattice regime |6| in the sense that there are one localized electron and one 
conduction electron per a unit cell. The ground state is described by a collection of local singlets 
between them. The efi^ective exchange coupling between the p- and d-sites J comparable 
to the hopping amplitude between p-sites. 

For I All ^ 1, A2 ^ 1, the correlations are suppressed. This is due to that the system decouples 
to a collection of the pairs of p- and d-sites. 

The ground state ( [4.57| ) is a half-filling state. The filling factor corresponds to that of a band 
insulator in the non-interacting system where the excitation gap satisfies the relation 1 < A < 2. 
Therefore, the correlation length is finite and is almost independent on the parameters Ai and A2, 
which is different from that of the ground state ( 4.57] ). The properties of the ground state (4.57) 



are completely different from that of the non-interacting system. 

D. Model C 

1. Hamiltonian 

The lattice of Model C is constructed from a cell with four d-sites and one p-site (Fig. |35|(a)). 
The model has four free parameters. We investigate the simplest model with one parameter. The 
cell Hamiltonian is obtained by choosing al^} = J2r=i Kcr,(7 = Ai cf^^ +A2 C2^^ +A3 C3 +A4 cf^^ 
+A5 C5 ^ and setting Ar = 1 (r = 1, 2, 3, 4) and A5 = A (see Fig. |3^(a) for intra-cell index). The full 
Hamiltonian is obtained by identifying sites 1 and 2 in the (n — l)-th cell with sites 3 and 4 in the 
n-th cell, respectively (Fig. |35|(b)). The Hamiltonian is 

N 

H^ = V V I V [ f -Ac^'^+cP _ Ac-^^t p _ ;^ dif p _ x d2^ p 

-"•S ' I [ \^ ^ ^71, o- '-71,0- ■^'^n,a'^n,a ^'^n+l,a^n,a ■^'^n+l,a'^n,a 

cr=t,i n=l 

_ dlf d2 _ rflf dl _ d2] d2 _ dl\ d2 _ d2\ dl , 
"-n,t7"-'n,tT "-n+l,cr"'n,(T ^n+l,a^n,a ^n-\-l,a^n,a ^n+l,a^n,cr ' "<.<-. 

I frPpPf pP t ^d dl] dl d d2] d2 

+ ^AT+l'^^+l.crCiV+l.o- + £Ar+lC^+l,crC^+l,o-}^> (4.78) 

where the on-site potentials are = — A^, = —4 (2 < n < A^), and ef = eJv+i = ~2. A unit cell 
is labeled by n. Here (cf^a) is the annihilation operator on a d-site for r = 3 (4) in the n-th 
cell. The ground state is 



9F, 



N 



C 



^3.1=^11 n « 



n,<T 



N 

-p TT TT f ,„d2^ , dij d2 1 , \ „p t 



, 



(4.79) 



which is a l/S-filhng state. 



2. Band structure in the single- electron problem 



We investigate the single-electron problem for the Hamiltonian ( [4.78 ). We consider the system 
with an even number of the cells under the periodic boundary conditions. The similar calculation 
to that in Sec. [V A 2| leads the dispersion relations 



El = 


1 " 

~2 . 


E2 = 


-2 


E3 = 


1 " 

~2 



6 + + 4 cos A; + J (6 + A2 + 4 cos ky - 8A2 



6 + A^ + 4cosA;- J (6 + A^ + 4cos A;) - 8A2 



(4.80) 



where 1, 2, and 3 are the band index, and k is the wave vector with ( [4.6[ ). The energy gap between 
the lowest two bands is A = for < A < a/2 and A = A^ for A > a/2. 

The electron number in the ground state (|4.7S| ) corresponds to full- filling of the lowest band. 
Therefore, the ground state of the non-interacting system is metallic for < A < -v/2 and is 
insulating for A > \/2. 



3. Correlation functions 
Prom (|]5|), the norm of the ground state is 



c 

G.S. 



^G.s.) — ~ {RiiLn + AR21L21 + -Rsi-^^si) 



ci 



where the corresponding matrices in ( [4.60 ) are 



D 



/ 2 + 4A2 + A^ 4A2 + 4A4 A^ \ / 1 

2 + A2 2 + 4A2 a2 , /= 

2 8 2 / \ 

/eiOO\ /ciOO 

62 I , and C 

V 63 



F 



(I 

4 



C2 
V C3 



(4.81) 



(4.82) 



The matrix T„ and the initial and the final vectors are derived in Appendix HI. Here (z = 1, 2, 
and 3) are the eigenvalues of T„ 



1 / t 
ei = - s + - 
3 V q 

s i 

62 = 1 

3 12 



63 



I 

12 



(1 + 2^/3) ^ + (1 - 2\/3) q 
1 + 2^/3) ^ + (1 - 2\/3 



9fi 



where s = 6 + SA^ + A^ t = X^{A8 + 58X^ + 16A^ + A^), and 



pi + 32p2i ) with pi 



468A^ + 512A*5 + 183X^ + 2AX^° + A^^ and p2 = A^(4096 + 6736A2 + 4288A^ +1328A6 + 192A^ + 9X^^). 
They satisfy ei > 62 > 63 > for A 7^ 0. The matrix L = (Lij) (R = {Rij)) is constructed from 
the left(right) eigenvectors. We choose them as we did in Sec. |V^, where Lji = {ej - 2)^ - 4ej A^ 
Lj2 = AX'^iejX'^ + ej - 2), Ljs = X'^{ej + 2), Rij = {ej - 2)^ - AejX'^, R2j = (2 + X'^)ej - 4, and 
R3j = 2(ej + 6). 

We evaluate the expectation value of the number operator {nf^). The transfer matrices asso- 



ciated with nf^ are 



N 



2A2 + A4 2A2 + 4A4 A^ 
- I 1A2 2A2 iA2 





N 



(d) 



/ 1 + A2 2A2 + A^ iA^ 
i l + A^ iA^ 



(4.83) 



V 



2 1 

They are derived in Appendix^. From (^^, (|^, (|48^ ), and (^^8^) , we obtain 



(4.84) 



^(4Liii?2i + ^12^31) for a = p 

^ ^ [(2Lii - L12 + 2L13) + (2Li2 - 4Li3) R21 + 2Li3i?3i] for a = d. 

It can be verified that there are two electrons per a unit cell as (n^) + 2{nf) = 2. The results are 
shown in Fig. ^ For |A| <^ 1, on-site potentials satisfy the relations < ep and |ep — e^^l ^ |A| 
(the hybridization). There is almost one electron per d-site. The p-sites are almost empty. For 
A 3> 1, on-site potentials satisfy the relations ep <C and — Cp ^ A. The p-sites are almost 
doubly occupied. The d-sites are almost empty. 

From (|Ol|) , (|4.32|) , (g^) , (|]8l|), and (|4.83D , the density correlation functions are 



where 



«>«) 



i?f ) = ALnRij + Li2i?3i 
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iST^-^s ^3 loi^ " - P 
iHI^^s ^3 foi a - a, 
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ip) 



J 4Ljii?21 + Lj2Rzi 

R^p = (2Lii — L12 + 2L13) + 2 (L12 — 2L13) i?2j + 2Li3i?3j 

-^j'^^ = -^jl (2-Rll) + Lj2 { — Rii + 2ii2l) + 2Lj3 (i?ii — 2i?21 + R?ji) ■ 

The density correlation functions take negative value and decay exponentially with distance. We 
show the results, in Fig. 38. For |A| <C 1 and A ^ 1, the density correlations are suppressed. 
We evaluate the spin correlation function. The transfer matrices are 



(lAMA^, 0), S 



Md) 



\X^, 0, 



They are derived in Appendix H3. From (^4.6 

A4 



lL,(d) 




(||8TD, and ( |486|) , we have 
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ei 



2^ {LiiRii + 5L12-R21) for /3 = p and a = p 



A2 



LiiRii 



'2ib7^i3^ii 



for (3 = p and a = d 
for f3 = d and a = d. 



(4.86) 



(4.87) 
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The spin correlation functions decay exponentially with the oscillating sign. We show the results 
in Fig. For |A| <C 1 and A » 1, the spin correlations are suppressed. 

We evaluate the singlet-pair correlation function. For ( [4.55 )-(iv), the correlation functions are 
vanishing for |/c — i| > 2. We evaluate them for (f4.55| )-(i) and (iii). (The position of the operators 
are shown in Figs, ^(e) and (f) in Appendix H4.) The transfer matrices are 

/-4A2 



2, H 



A' ( 1, 5, ) 



H 



V 








(4i 



for ^ = p and dp. They are derived in Appendix H4. From ( [4.68| ), ( 4.81 ), and ( 4.8^ ), we obtain 



eici 



(Liii?ii + 5Li2i22i) for a, (3 = p, dp. 



(4.89) 



The singlet-pair correlation functions decay exponentially with distance. We show the results in 



Fig. 40. For |A| <C 1 and A ^> 1, the correlations are suppressed 

ransfi 

1 + A2 A2 



We evaluate the correlation function (q^Cj^). The transfer matrices are 



-A 



2 2 
1 4 1 



G 



R,{d) 



1 + A2 A2 
1 1 



1 



QL,id) 



1 




(4.90) 



A^ A^ 
i + A2 A2 



They are derived in Appendix H5. The matrix G„ is diagonalized as 

G„ = -2RDC^L, 



(4.91) 



where the matrices are shown in ( 4.71| ). The corresponding quantities are gi = (2 -|- A^ -|- Xuj^)/2 
and 52 = (2 - -f-^Aa;3)/2 with ^3 = \/A + A^. They satisfy gi > g2> for A / 0. We choose 
the matrices L = {Lij) and R = {Rij) as we did in Sec. IV where Lji = {gj — 1)/A^, Lj2 = 1, 
Rij = gj — 1, and i?2j = 1- (See also (4.72).) From these matrices and (4.44), we obtain 



where 



[2(5. 



-2^ 
ei 



ft' + 



-2^ 

ei 



(a) 



(4.92) 
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'^ciCmeig 



l)-^^ll + gniL 12] 

X 2{g^ - l){Rn + R21) + X^Rii + 4i?2i + i?3i) 



id) 



l)Li2 + 4L 



13j 



{gm - l)i?ll + X^R21 



(4.93) 
(4.94) 



The correlation functions decay exponentially with the oscillating sign. 

By the Fourier transformation of the correlation function (Cj^cj^), we obtain the momentum 
distribution functions for a = p and d 

(4.95) 



K.) = /^Vi(A:,ri) + ft'F2{k,r2) + ft\ 



where Fi{k, Vi) is defined by ( [4.77 ), rj = —gi/ei, and /g"^ = (nf^). The results are shown in Fig. |4l| . 
There is no singularity in the momentum distribution functions. For the momentum distribution 
for j)-site, it has a sharp peak at /c = tt for |A| <C 1. The peak is, however, expected to vanish in 
the complete limit. It is almost fiat for A S> 1. The momentum distribution for d-site is expected 
to be flat for the complete limit |A| <C 1. 
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4- Discussion 



All the correlation functions under consideration decay exponentially with distance. These 
results suggest the existence of a finite excitation gap. Therefore, it is expected that the state 



is insulating. 



In 



Their correlation lengths are ^,1 
1 



In 



6 



(,bb, = 



In 



and 



(Fig. |4^). (They satisfy the relation ^cc > ^nn > S,ss = ^bb-) We note that the 
spin correlation is antiferromagnetic. 

We consider the region |A| <^ 1. There is almost one electron per a d-site. The p-sites are 
almost empty. The correlation lengths are large. However, the correlation functions for the nearest 
neighbor sites vanish. 

For A ^ 1, There is almost two electrons per a p-site. The d-sites are almost empty. The 
correlation lengths and the correlation functions for the nearest neighbor sites vanish. 

The ground state ( [4.79 ) is a 1/3-filling state. The filling factor corresponds to that of the band 
insulator in the non-interacting system for A > \/2. In the non-interacting system, we have the 
metal-insulator transition at A = \/2 by the variation of A. However, the ground state ( [4.79| ) is 
insulating for any A. The properties of these states are completely different. 



9Q 



V. ABSENCE OF THE PERSISTENT CURRENT 



In Sec. [V, we have assumed the parameter to be real. Relaxing the condition, effects of a 
magnetic field is included by taking hopping matrix elements to be complex. Thus the effects of 
the magnetic field is investigated exactly for the systems of strongly correlated electrons described 
in this paper. Let us calculate the persistent current. Considering the system in a ring geometry 
and putting a flux through the ring, we can measure the Aharonov-Bohm effect and the persistent 
current ,2^,21,22,23|. In the ring geometry, we include the effect of the flux $ by changing the 
hopping matrix elements of the A^-th cell. We first classify sites in the A^-th cell into two classes as: 
(i) sites which belong to the A^-th unit cell and (ii) sites which are identified with sites in the cell 
Ci. We denote the sets of the sites (i) and (ii) by Cn-u and CN-e-, respectively. The cell Hamiltonian 
(P^) associated with the iV-th cell is obtained by choosing 



\Cn\ 



where 



r=l 



A^^^e^* for r E CN-e 



with real A 



{N) 



. From (|0|) and 



^ Ar^^ for r G Cn-u, 
for 1 < n < — 1, the Hamiltonian is 

<T=^,l x,j/eAjv 



where 



N-l 



n=l 



($), for X = f{n, r) and y = f{n, s), 



with 



^r,s 



AN) 



($) 



At As On 

2(a1") " 
f (Ar)(cl>) 

^{N) T(iV), 
Ar As \ 

Ar As 

2(a1^)^' 

a( 



r 

(a(^)) 



for r ^ s 

for r = s and r G Cn-u=oo 
for r = s and r G C„;[/=o, 

(l<m<Af-l, l<n<N 
for r 7^ s and r G Cn-c-, s G Cj\f;u 
^) for r ^ s and r G Cn;u, s G CN;e 
for r 7^ s and r, s G Cj\f;u 
for r = s and r G Cn-u=oo 
for r = s and r G Cn-u=o- 



1) 



AN) 



The ground state of Hs{^) is 



G.S. 



V n 



n 



a. 



a 



N. 



The ground state energy is given explicitly by 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



N-l 

n=l rGCn tGCn 
N 

= -EE2|a(") 

n=l reCn 



(5.8) 



It is independent of the flux <I>. The persistent current / is evaluated by using the Byers-Yang 
relation [lyl. We obtain 



/oc 



0. 



(5.9) 



The persistent current is vanishing for any of the solvable models discussed here. This is consistent 
with our conclusion that the ground state is insulating. 

Extending the discussion here, the absence of the persistent current can be shown in any 
dimensions for the models discussed in this paper. 



VI. SUMMARY 

We investigated three models with strongly correlated electrons which have the RVB state as an 
exact ground state. The number of the electrons per unit cell is restricted to be 2. The correlation 
functions are evaluated exactly using the transfer matrix method for the geometric representations 
of the valence-bond states |^^. The two-point correlation functions for spin, density, and singlet- 
Cooper-pairs are obtained for any distance. All the correlation functions decay exponentially with 
distance. The momentum distribution functions are also evaluated and there is no singularity. The 
results suggest that the ground states of the models are insulating. The persistent currents are also 
considered and turned out to be vanishing. 
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APPENDIX A 



We describe the models by several authors ||T|,0^^ and those investigated in this paper using 
the cell construction of Tasaki (see Table ^. The lattice is constructed by the cell in the second 
column in Table |. The cell for the line graph is constructed as follows. We define a lattice 
£ = (A, B) where A is the set of the sites (vertices) and B is the set of the bonds (edges). The line 
graph L{C) = {A^,B^) constructed from a lattice A has the bonds of C as sites (A^ = B), and 
two sites are connected by a bond in B^ if the corresponding bonds in B have a site in common. 
The cell is defined by sites (g A^) which are connected by a bond (g B) to a same site (e A). 



APPENDIX B 



We first note the equalities 



i^a^lj = sgnia) c]^^ [via) 



and 



(Bl) 
(B2) 



where Aj = -Y.a='l,l4,o- (4-(7S','^) ^i-a which has the property Ai|0) = (0|Ai = 0. 

We show an equality. Consider a connected graph W U W with 2n (n > 1) bonds. We set 
W = {{2, 3}, {4, 5}, • • • , {2n, 2n + 1}} and W = {{1, 2}, {3, 4}, • • • , {2n - 1, 2n}}. We define the 
quantity 



L{n;a,p) = ( 



2k 



, k=l 



k=l 







For n = 1, from (|B1| ) we have 

L{l;a,p) = (O 



b2A,P (^1X2) I O) 



-sgn{a) (o 



^2,34,p'^2,-cr 



sgn{a)sgn{-p) (o 



^2~a^2-a^Z,p^?,,p 



Prom (iBll), we have 



L{n;a,p) = ( 



II ^2fc,2fc+l ) '^l,(74n+l,p4,2 H 4fc-l,2fc 



, k=l 

-sgn{a) ^0 
-sgn{a) ( 



k=2 



II ^2fc,2fc+l )4n+l,p4,-cr H 4a:-1; 



2A: 



, A;=l 
n 



fc=2 



II ^2fc,2fe+l ]h,3'^2n+l,p'^2-a II ^2fc-l 



2fc 



, k=2 



k=2 



-sgn{a) sgn{—a) ( 



II ^2fc,2A:+l ( C3^o-)4n+l,p H 4fc-l, 



2fc 



. k=2 



k=2 



II ^2fc,2A:+l I '^3,(T4n+l,p H 4fc-l; 



2fc 



. A:=2 



k=2 



(B3) 



(B4) 



(B5) 



We change the labehng of the lattice sites by the rule j ^ j — 2 and obtain 

'n-l \ 



L(n;cr, p) 







\ n-l 
n ^2fc,2fc+l ) '^l,o-'^2n-l,p n ^2fc-l,2fc 



. k=l 

-L{n- l;a,p) 



k=i 







(B6) 



We show ( ^j.lOl) . For the self-closed bonds and degenerate graph we have Wj = 1 and Wj = 2, 
respectively. Consider the non-degenerate loop Wj U Wj with Ij = 2n (n > 1) bonds. We set Wj = 
{{2, 3}, {4, 5}, • • • , {2n - 2, 2n - 1}, {2n, 2n + 1(= 1)}} and Wj = {{1, 2}, {3, 4}, • • • , {2n - 1, 2n}}. 
From (|B2| ) and ( |B6| ) the weight is 



II ^2fc,2fc+l II ^\k-l,2k 
k=l k=l 
/n-l \ 

^ II ^2fc,2A;+l J ^2n,2n+ 
/n-l \ / 



n-l 

1^2n-l,2n II ^2k-l,2k 
k=l 



n-l 



n ^2fc,2fc+l J2 C2n,a4n,aCl,-a4n-l,-a + ^Sn H 4fc-L 



2fc 



, fc=l 



k=l 



E 



'n-l 



n-l 



II ^2fc,2fc+l I '^2n,(74n,(7Ci _o-C2n-l,-cr II 4fc-l,2fc 



. k=l 

•n-l 



fe=l 



n-l 



~ X! ( II ^2fe,2fc+l ) "^1 -(74(n-l)+l,-(T II 4a:-1,2A: 

cr=1,l \ \ k=l ) k=l 

= 2(-l)"-^ 
= 2(-l)'^/2-i, 

where we used the relation / = 2?7- in the last line. 



(B7) 



APPENDIX C 



We first note the equalities 



where 



Sfk 



Sfb. 



4i - 4,T4a + 4a4,T' 



is the creation operator of the triplet-pair between sites i and j and satisfies the relation b. 
Here 



(CI) 



(C2) 



kj = 








2 


kj = 






2 



(C3) 

which have the properties Si\0) = {0\5i = and 6i\0) = {0\6i = 0. 

We show (|3.17] ). Consider the connected graph W^^'y'^ U W'^^'^^ with 2n (n > 1) bonds. We 
set = {{2, 3}, {4, 5}, • • • , {2n, 2n + 1}} and W'^'^'^'^ = {{1, 2}, {3, 4}, • • • , {2n - 1, 2n}}. We 

denote the right hand side of ( p-lT] ) by S{m) as 



S{m) = ( 



, fc=l 



2A; 



fc=l 



, 



(C4) 



where we set x = 1 and y = m. From (pi|), we have 



S{2m + 1) = ( 



vfc = l / fc = l 

' m— 1 

n ^2fc,2fc+l 
V fc=l 







m—1 \ 

n ^2fc,2fc+l ) 
fc=l / 

m—1 \ 

n ^2fc,2fc+l ) 
fc=l / 

/ m—1 

n ^2fc,2fc+l 



fc=l 
n 



n ^2fc,2fc+l 
lc=l 
n 

n ^2fc,2fc+l 
k=l 

n 

II ^2fe,2fe+l 
k=l 

n 

II ^2fc,2fc+l 



, fc=l 



II ^2fc,2fc+l j '^2m,2m+ 
A:=m+1 / 

1' 



I'S'lm+l'S'l H ^2k-l,2k 
k=l 







II ^2fc,2fc+l ( 9^2m+l,2m + 



''2m+l,2m 



^1 II 4fc-l,2fc 



fc=m+l 
n 

II ^2fc,2fc+l ) ( ~9^2m,2m+l ) ^1 II ^2fc-l 
fe=m+l / 

/ n 

\ II ^2fc,2fc+l 
\ A;=m+1 
n 



k=l 





2k 



k=l 
n 



^2m-l,2m^2mSl H ^2fc-l,2fc 



fc=l 



'S'l'S'lm H ^2fc-l,2fc 



fe=l 







' m—1 



'S'l'S'2m^2m-l,2m I II ^2fc-l,2fc ) II ^2fc-l,2fe 



fc=l 



A:=m+1 
' m—1 



•S"! f r,^2m,2m-l + '^2m,2m-l ) [ H ^L-l,2fc ) II ^L-l. 
^ ^ \ k=l / fc=m+l 

"S"! (~2^2m-l,2m) ( II 4a:-1,2A: ) II 4fc-l,; 
^ ^ \ k=l / fc=m+l 



2A; 



(-1)^ 



(-1)^ 



' m—1 



II ^2fc,2fc+l '^l'^2m-l4m-l,2m II 4fc-l,2fc II 4fc-l, 



2k 



. k=l 



. k=l 



k=l 



k=m+l 



2fc,2fc+l I'S'l'S'lm-l H ^2k-l,2k 
/ k=l 



{-fS{2m-l) 



{-ly^Sil) 

(_l)2m^0 



II hk,2k+l jSfSf Y[ 4fc-l, 



2k 



. k=l 

n 



k=l 



II ^2fc,2fc+l II ^2fc-l,2fc 



k=l 



k=l 



(C5) 



.^4 



Similarly, we have 



5(2m) = (-l)2'"-i^(0 



k=l k=l 



, 



Therefore, we obtain 



5(m) = ^{-l)'^^'''y^Wj 



where we used the relation d{x, y) = m — 1 and the definition of Wj. 

APPENDIX D 

We show (13321) . Consider the line U W'^'^''^'^ with 2n (n > 1) bonds. We set 

{{2, 3}, {4, 5}, • • • , {2n, 2n + 1(= y')}} and TV' ("''^^'^ = {{!(= x'), 2}, {3, 4}, • • • , {2n - 1, 2n}} 
), the weight is 



n ^2k,2k+l ) '^l,o-'^2n+l,p II ^2fc-l,2fc 
, fc=l / A:=l 



Oy = L(n; u, cr) 

= (-ir 

= (_i)'(^'y)/2 



where we used the relation l(x',y') = 2n. 



APPENDIX E 



We show ( p8|) . Consider the graph W^'^^UW'^^^ we set M^(^) = {{l(= x), 2}, {3, 4}, 
1, 2n}} and W'^""^ = {{2, 3}, {4, 5}, • • • , {2n, 1}}. We have 



n ^2k-l,2k I "-1,(7 n 4fc,2fe+l 



, k=l 



k=l 



n— 1 



n hk-l,2k )^1,2 (l C]^,o-4,cr) 4n,l II 4^,2^+1 



. k=2 



k=l 



II ^2k~l,2k II ^2k,2k+l 
k=l k=l 

From (|B3|), the second term is 



0-0 



II ^2k-l,2k ] "^1,0-^1,0- H ^2fc,2fc+l 
, fc=l / fc=l 



. 



, k=l 



II ^2fc-l,2fc ) '^1,0-^1,0- H ^2k,2k+l 
/ A:=l 
n \ 7 

II ^2k-l,2k ]^l,2'^l,a'^l,cr^2n,l H ^2fc,2fe+l 
fc=2 / fc=l 

n \ n— 1 

n ^2fc-l,2fc (-^l,lC2,a)(-4n,a4,l) 11 4fc,2fc+l 
t=2 / 

n \ n— 1 

II ^2k~l,2k ]^2,u^2n,a W ^2fe,2fe+l 



A:=l 



, k=2 
■n-l 



k=l 



\ n-l 

II ^2fc,2fc+l j '^l,cr^2n-l,cr II ^2fc-l,2fc 
A:=l / fc=l 



L{n — 1; (J, fj) 
^_-^y(x)/2-i_ 



(E2) 



The fourth hne is obtained by renaming the site j to j — 1. We used the relation l{x) = 2{n — 1). 
From (|B^) and (p]), we have 







n ^2fc-i,2fc^i,<7 n 



2fc,2A,-+l 



k=l 



k=l 



Oj = 2(_l)'W/2-l _ (^_iyix)/2-l 

= (_iyw/2-i_ 



(E3) 



APPENDIX F 

We show (13.351) . For (i), the sets M/^^'^^) and ly'^^'^^ are empty. The right hand side of (|^) is 



Wj{x,y;cr) = ^ {(o 



0) + (0 



1, 



(Fl) 



where we used the equality (0|c^_^Cj^_^4^^4^^|0) = {0\c^^^Cy_^cl_^cl^JO) = 1. 

For (ii), we suppose that the site x belongs to the graph. From (|3.28|) , the right hand side of 
O is 



Wj{x,y;a) 







^y, — a^y^ — a 



;he equality ( 
For (iii), the right hand side of (|3.34D 



{u',v'}ew 

o\/o 



( n j^u',v''^x,a(i:,a H ^u,v 



(F2) 



where we used the equality (O|cj^^o-cj JO) = {0\cy^_^Cy_^\0) = 1. 



IS 



n 



^u'v' 1 ^x,a^x,a 



' n ^i,. 



+ 



( n ^'y)vL n Co) 



(F3) 

where we used (|3.2^ ). 

For (iv), consider the graph VF^^'?') U W'^^'^y We set VF^^'f) = {{1(= x), 2}, {3, 4}, • • • , {2n - 
1, 2n}} and W'^^^'^^ = {{2, 3}, {4, 5}, • • • , {2n, 1}}. We denote one of the terms in the right hand 
side of ( |3.34 ) by G{n, j; a,T) as 



G(n,j>,r) = {0 
From (|BlJ ), we have 



. fc=l 



G(n,j;(T,r) = (o 



fc=i 



n-l 

n ^2fc-l,2fc \ ^l,2^l,a^j,T'^j,T^l,a^2n,l W_ ^2fc,2fc+l 
t=2 / fc=l 
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t=2 / 

71 \ n— 1 

n ^2fc-l,2fc ]^2,a'^j,T^j,T^2n,a W_ ^2k,2k+ 
k=2 / fc=l 

n— 1 \ n— 1 

n ^2fc,2fc+l ) '^l,o-S-l,-r'^j-l,r'^2n-l,cr II ^2fc-l, 
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(F4) 



n-l 



A;=l 



, fc=l 



fc=l 



(F5) 



The last hne is obtained by renaming of the lattice sites by the rule k ^ k — 1. For j = 2m, we 
have 



G{n - 1,2m; a, t) = ( 



■n-l 



n-l 



n ^2A:,2fc+l ]^l,cr^2m-l,T'^2m-l,T^2(n-l)+l,a 11 ^2k-l 
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We change the labeling of the lattice sites by the rule j — > j — 2 and obtain 

■n-2 \ n-2 

II ^2k,2k+l )'^l,o-C2(„,,_i)_i,r'^2(m-l)-l,T'^2(n-2)+l,(T II ^2k-l,2k 
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G{n - l,2m;a,T) = - ( 
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II ^2A:,2A:+1 ) '^l,(T'^l,r4(n-m-l)+l,(T ( 4,r4,2) II 4a:-L 



2k 



k=l 



C n—m—l \ 
n ^2fc,2fc+l j Vl,A,lC. 



-l)™sgn(T)5,,_, 
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^2fc,2fe+l 1^3 _r'^2(n-m-l)+l 
k=2 / k=2 
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?57 



We change the labehng of the lattice sites by the rule k ^ k — 2 and obtain 

■ n— m— 2 \ n— m— 2 



G{n - 1, 2m; a, r) = (-1)'"(5^,_^ ( 



2k,2k+l I ''l,-r''2(n-m-2)H 



n b., 

k=l 

= (-l)™<5^,_^L(n -m-2;a, -r) 

= X (-1)"— 

where we used ( |B6| ) . For j = 2m — 1 , the similar calculation leads 

G(2m-l;a,T) = (-l)"-i<5,_,. 
Combining @ and (||) is 



1,. n 

fc=i 



-l,2fc 







n ^2fc,2fc+l ]'^l,a^y,T'^y,T'A,a II 



2fc 



, fc=l 

We obtain the weight 



k=l 



=(-ir-'5., 



(_l)<i{a:,j;)^. 



(a;, y;cr) = ^ [G{n, y; a, a) + G{n, y; a, -a)] 

(^_iyix,y)/2-l 



(F8) 



(F9) 



(FIO) 



(Fll) 



where we used the relation n = l{x,y)/2. 
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APPENDIX G: DERIVATION OF TRANSFER MATRICES FOR MODEL B 



In Model B, the sites which is identified in the cell construction are p-sites and there exists the 
sites with four bonds in the geometric representation. We need the procedures shown in Fig. ^ (b) 
and (c) to calculate the contribution from the graph. Since we do not need the procedure shown 
in Fig. m (a) except for the correlation function (q.^c^^), mj = in ( 3.11[ ), ( p. 19 ), ( 3.30D , ( 3.37 ), 
and ( p^ ). 



1. Norm of the ground state 



For the sake of convenience, we draw the lattice shown in Fig. 25(b) as Fig. ^(a). The ground 
state is written by ( |3.3| ) geometrically, where an example of the valence-bond configuration V is 
shown in Fig. ^(b). The geometric representation of the norm is ( p.ll| ), where an example of the 
graph VuV is shown in Fig. |4^(c). We first evaluate the contribution from a graph VuV. It can 
be decomposed into the subgraphs UiUUl (z = 1, 2, • • •, niVUV')). The Loops which are extending 
over more than two cells are allowed. (See Figs, ^(c) and (j)). After the elimination of sites with 
four non-closed bonds (by the procedure shown in Fig. |5|(b) or (c)) we have one degenerate loop 
and self-closed bonds (Fig. |4^(d)). Therefore, the loop (Fig. |3| (j)) has weight 2. 

We derive the matrix T„ in ( [4.82| ). A subgraph Ui U is constructed by the five kinds of cells 
shown in (Fig. |4^ ). The sum over the graph VUV' is equivalent to that over all the combination of 
above five cells under the restriction that a p-site has at most four valence-bonds. (The restriction 
means, for example, that the identification of the right p-site in Fig. ^(b) with the left p-site in 
Fig. 0(c) is forbidden.) We shall always take into account the restriction hereafter and it should 
be understood implicitly. We have to distinguish three cases due to the restriction. Let j4„, Bn, 
and Cn be the quantity defined by the right-hand side of (|3.6D on the lattice A„. For An, the sum 
is taken over all the combination of the cells shown in Fig. ^ with the restriction that the n-th 
cell is represented by Fig. |4^(a) or (b). For Bn and C„, the sum is taken as was done for An with 
the restriction that the n-th cell is represented either by Fig. p^(c) or (d) for Bn and by (e) for Cn- 
They are represented diagrammatically 

An = (Gl) 

Bn = (G2) 
Cn = (G3) 

The recursion relations are 

A-n = 

= XiAn-i + 2Ai^„_i + 2XlBn-i (G4) 

Bn = 

= XfAn-^l + 2XlBn-l + XfCn-l + XfXlAn-l + A^A^-Bn-l (G5) 
Cn = 

= XfAn-1 + 2X'lBn-l + XfCn~l, (G6) 

where we adopt the following rule to assign weight 2 to the loop which is extending over more than 
two cells. We assign the coefficient 1 for all the terms in Bn, since the loop may continue to the n+1- 
th cell. In the third term in An and the second term in Bn and C„ we assign weight 2, since we make 
sure that the loop finishes there. The corresponding matrix is T„ in ( 4.601 ). Since any cell in Fig. 



44 is allowed at the boundaries, the initial and the final vectors are / = {Aq, Bq,Co)'^ = (1,0,0)^ 
and F = {An,Bn, Cn)'^ = (1, 2, 1)'^, respectively. 



2. Expectation value of the number operator 



The geometric representation of the expectation value is ( |3.29| ), where an example of the graph 
V U V is shown in Fig. ^(e). We derive the matrix n|^\ We have five kinds of graphs on the 
i-th. cell (Fig. |45|(a)-(e)). Let a[^\ and Q-^'' be the quantity defined by the right-hand side 

of ( 3.29| ) on the lattice Aj. The restriction for the sum is that the i-th cell is represented either by 
Fig. ^(a) or (b) for A^l'\ by (c) or (d) for and by (e) for C^^\ The recursion relations are 

i(p) _ 



X^A.^i + 2A^5,_i + A^Ai_i (G7) 



-A^Ag^i-i + A^A^-Bj-i + A^i?j_i + \\Ci-i (G8) 



\\B^-l, (G9) 



where we used the same rule in subsection p 1| for the subgraph without the number operator. 
For the subgraph with the operator, we adopt the following rule to assign the weight ( p. 28 ). For 
the first and the third terms in the third and fourth terms in and c\^\ we assigned 

coefficient 1. For the second term in A^^\ we assigned weight 2, since the loop is decomposed 
into self-closed bonds and a degenerate loop which finishes there. For the second term in Bl \ we 
assigned coefficient 1, since the loop is decomposed into self-closed bonds and a degenerate loop 
which may continue to the i + 1-th cell. For the first term in B'f \ we assinged coefficient 1/2 to 
cancel weight 2 which is assigned at another end of the loop where we regard it as the end of a loop 
without the number operator. The corresponding matrix is N-^^ in ( 4.61| ). Using the cells shown 
in Fig. E^(f) and (g), the similar calculation leads N-'^'' in (4.61). 



3. Spin correlation function 

The geometric representation of the expectation value is ( |3.18| ), where an example of the graph 
V UV is shown in Fig. ^(f). From the procedure in Fig. |5|(b) or (c), the graph with the spin 
operators is decomposed into a degenerate loop with the sites i and j and the self-closed bonds. 
Therefore, d{i,j) = 1 in ( 3.1S| ). The degenerate loop has weight —1x^x2. We assign ^ at the 



i-th cell and — ^ x 2 at the j-th. cell in the recursion relation. We derive the matrix S„ in (4.54) 



Let Sn be the quantity defined by a sum. The sum is taken over VuV on the lattice A„ such that 
the graph consists of loops shown in Fig. ^ on the fe-th cell (l<A;<i — 1) and that shown in Fig. 
|4^(c) on the l-th cell [i <l <n). The recursion relation is 

Sn — 

= \l\lSn-l, (GIO) 



and the corresponding matrix is S„ in ( 4.64 ). 



We derive the matrix S^'^^^ The graph on the i-th cell is shown in Fig. H6|(a). Let be the 



quantity defined by a sum. The sum is taken over V UV on the lattice Aj such that the graph 



consists of loops shown in Fig. 44 on the A;-th cell (1 < A; < i — 1) and that shown in Fig. |46|(a) on 
the i-th cell. We have 



C(P) 



^ X XlxlA^i. (Gil) 



4n 



The corresponding matrix is S^'*^^^ in (4.64). 



We derive the matrix S 



L,{P) 



c. 



ip) 



Let 5]^^ and 



J be the quantity defined by the right-hand side of ( |3.1^ 
for the sum is that the j-th cell is represented either by Fig. 
recursion relations are 



The graph on the j-th cell is shown in Fig. p6|(b) 

m the lattice Aj. The restriction 
|(d) for Bj and by (e) for Cj. The 







X 2 X XlSj^i 



(G12) 



(G13) 



(p) 



1 



X 2 X \\Sj-i. 



The corresponding matrix is S^'^^^ in ( [4.641 ). From the graph shown in Fig. 



(G14) 

1(c) ((d)), the similar 



calculation leads the matrix S„- 



4. Singlet-pair correlation function 



The geometric representation of the expectation value is (3.11), where an example of the graph 
y U y is shown in Figs, ^(g) and (h) for ( [4.55| )-(i) and (iii), respectively. In the Figures, a double 
solid (broken) line represent the operator b\ ■ (bk,i)- We derive the transfer matrix H„ in (4.54). 
Let Hn be the quantity defined by a sum. The sum is taken over V U V on the lattice A„ such 
that the graph consists of loops shown in Fig. ^ on the fe-th cell {1 < k < i — 1) and graph shown 
in Fig. ^(a) on the l-th cell {i <l <n). The recursion relation is 

Hn = 



A2H^ 



(G15) 



and the corresponding matrix is H„ in ( [4.67] ). 



We derive the matrix H^'^^^\ (We consider the case (4.55)-(i).) The graph on the i-th cell is 
shown in Fig. ^(b). Let h\^^^ be the quantity defined by a sum. The sum is taken over V [JV on 
the lattice Aj such that the graph consists of loops shown in Fig. 44 on the n-th cell (1 < n < i — 
and that shown in Fig. 47(b) on the i-th cell. The recursion relation is 



(pp) 



1- 



(G16) 



The corresponding matrix is H 



R,(pp) 



m 



We derive the matrix H 



L,{pp) 



(We consider the case (|4.55D -(i).) We have two kinds of graphs 



on the fc-th cell (Fig. 47(c) and (d)). Let B^^^^ and Cf^^' be the quantity defined by the right-hand 
side of ( |3.40| ) on the lattice . The sum is taken over V UV with the restriction that the k-th. 



(pp) 



cell is represented by Fig. ^(c) for B^f^^ and (d) for C^^\ The recursion relations are 



(pp) 



A 
B 



(pp) 

k 

(pp) 



(pp) 



2Ai//fc_i, 



(G17) 



(G18) 



(G19) 
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The corresponding matrix is Hjf^) in ( pTD - 

We consider the case ( 4.55| )-(iii). The graphs on the i-th and the k-th cells are shown in Fig. 
(e) and (f), respectively. The similar calculation leads the matrices H^'^'^^^ and H^''-'^^^ in (4.67). 



5. Correlation function (Cj^c^^) 



The geometric representation of the expectation value is ( | 3.23| ), where an example of the graph 
VLIV is shown in Fig. ^(i). We derive the matrix G„ in ( 4.54|) . A line is constructed by combining 
three kinds of cells shown in Figs. p8|(a)-(c). To calculate the weight of the graph, it is necessary 
to use the procedure shown in Fig. ^(a). If we use the procedure m times for the graph with 2n 
bonds, we obtain a line with 2n — 2m bonds and m self-closed bonds. The weight of the line is 
(_l)(2n-2r?i)/2 ^ where the former power comes from ( 3.22 ) and the latter one comes from 

the procedure shown in Fig. |5|(a). We use the same rule as that in ( 4.37] ) and the minus sign in Fig. 
|^(a) is automatically taken into account. Let Dn and En be the quantity defined by a sum. The 
sum is taken over V U V' on the lattice A„ such that the graph consists of loops shown in Fig. ^ 
on the fe-th cell {1 < k < i — 1) and the graphs shown in Fig. |4^(a)-(c) on the l-th cell {i < I < n). 
The restriction for the sum is that the n-th cell is represented by (a) or (b) for D„ and (c) for En- 
The recursion relations are 



En 



-\l\2Dn-l — \l\2En-l 



\l\\Dn-l 



(G20) 



-\l\2Dn-l — \i\2En-l- 



(G21) 



The corresponding matrix is G„ in ( [4.70 ) 



We derive the matrix G^'''^^ We have three kinds of graphs on the i-th cell (Fig. |48|(d)-(f)). 



Let o'f^ and Ef be the quantity defined by a sum. The sum is taken over V L}V' on the lattice 



Aj such that the graph consists of loops shown in Fig. 44 on the /c-th cell {1 <k <i—l) and those 



shown in Fig. 48(d)-(f) on i-th cell. The restriction for the sum is that the i-ih cell is represented 
either by Fig. || (d) or (e) for and by (f) for E'f\ The recursion relations are 



(p) 



E. 



(p) 



The corresponding matrix is Gj in ( 4.70| ). 



-\\\2Ai 

R,ip) 



1- 



(G22) 
(G23) 



We derive the matrix G^'^^\ We have four kinds of graphs on the j-th cell (Fig. ^(g)-(j)). 



Let A^^\ and cj^'' be the quantity defined by the right-hand side of ( 3.23| ) on the lattice Aj. 
The sum is taken over V UV such that the graph consists of loops shown in Fig. ^ on the A;-th 
(1 < < i — 1) cells and a line shown in Fig. ^(a)-(c) on the l-th (i < ^ < i — 1) cells. The 
restriction for the sum is that the j-th cell is represented by either by Fig. p8|(g) for by (h) or 
(i) for and (j) for 



A 



The recursion relations are 

(p) _ 



A^Z^i-i 



(G24) 
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B 



(p) 



C. 



ip) 



j(p) 



(G25) 



(G26) 



The coefficient 1/2 in the first term in is to cancel weight 2 which is assigned at the left end 
of the line where we regard it as the end of a loop. The corresponding matrix is G^'^^'* in ( 4.7C| ). 
We derive the matrix 0^'^"^^ We have two kinds of graphs on the i-th cell (Fig. ^(k) and (1)). 



Let D^"^^ and E^""' be the quantity defined by a sum. The sum is taken over V UV on the lattice 
Aj such that the graph consists of loops shown in Fig. 44 on the A;-th cell (l<A;<i— 1) and that 



shown in Fig. 
either by Fig. 



l(k) or (1) on z-th cell. The restriction for the sum is that the y-th. cell is represented 



8|(k) for Dy^^ and by (1) for Ey""' . The recursion relations are 



id) 



id) 



-AiA2^i_i 



E.. 



id) 



(G27) 
(G28) 



— —XfAi^i — 2Afi?j_i — XfCi-i. 

The corresponding matrix is G^'^'^^ in ( 4.7C| ). 

We derive the matrix G^'^'^\ We have two kinds of graphs on the j-th. cell (Fig. 48 (m) and 

(n)). Let A^'^^ and Bj^^ be the quantity defined by the right-hand side of ( 3.23| ) on the lattice Aj. 
The sum is taken over V UV' such that the graph consists of loops shown in Fig. 44 on the k-th 
(1 < k < i — 1) cells and a line shown in Fig. ^(a)-(c) on the Z-th {i < I < j — 1) cells. The 
restriction for the sum is that the j-th cell is represented either by Fig. ^(m) for A^'^^ and by (n) 
for Bj'^\ The recursion relations are 



A 



id) 



B 



(d) 



id) 



-XiD,^ 



--XiX2Dj-i - -XiX2Ej-i 



0. 



(G29) 

(G30) 
(G31) 



The coefficient 1/2 in Bj'^^ is to cancel weight 2 which is assigned at the end of the line where we 



regard it as the end of a loop. The corresponding matrix is G^'^'^^ in ( [4.70| ). 



APPENDIX H: DERIVATION OF TRANSFER MATRICES FOR MODEL C 



In Model C, the sites which is identified in the cell construction are d-sites and there is no site 
with four bonds in the geometric representation. Therefire, we do not need the procedures shown 
in Figs. |, rrij = in (|3Tl|), {KM, (KM), (EM), and (KM. 



1. Norm of the ground state 



The ground state is written by ( |3.3D geometrically, where an example of the valence-bond 
configuration V is shown in Fig. |4^(a). The geometric representation of the norm is ( |3.11 ), where 



an example of the graph V UV is shown in Fig. |4^(b). We first evaluate the contribution from a 
graph V U V . It can be decomposed into the subgraphs Ui U {i = 1, 2, ■ ■ n{V U V')). 



We derive the matrix T„ in (4.82). A subgraph Ui U [/{ is constructed by 27 kinds of cells 
shown in Fig. 5C . The sum over the graph V dV is equivalent to that over all the combination 
of above 27 kinds of cells under the restriction that a d-site has at most two valence-bonds. We 
shall always take into account the restriction hereafter and it should be understood implicitly. Due 
to the restriction, we have to distinguish four cases with respect to two d-sites at left end of the 
lattice A„: (i) they do not belong to valence bond; (ii) one of them belongs to two bonds; (iii) they 
each belong to two distinct bonds; (iv) they belong to a single degenerate loop. Let An-Dn be the 
quantity defined by the right-hand side of (^|^) on the lattice A„,. The sum is taken over V UV' 
such that the graph consists of the subgraphs shown in Fig. ^ The restriction for the sum is that 
the n-th cell is classified by the condition (i) for An, (ii) for (iii) for C„, and (iv) for Z)„. They 
are represented diagrammatically 

An = (HI) 
Bn = = (H2) 
Cn = = (H3) 
Dn = (H4) 

For the weight of a non-degenerate loop with 2n bonds, we assign the sign part (—1)""^ and the 
coefficient 2 separately. The loop belongs to n cells, since a cell has 2 bonds. We assign —1 to each 
n — 1 cells except for the cell in the right end of the loop. Accordingly, every recursion from C„_i 
gives a minus sign. At left end of the loop, we assign the weight 2. Therefore, when we attach 
Cn-i to the n-th cell, we adopt the following rule. If the n-th cell satisfy the condition (i) or (ii), 
we make sure that the loop finishes on the n-th cell. Therefore, we assign coefficient —1 x 2. In 
other cases, we assign coefficient —1. The recursion relations are 



■An 



Bn 



Cn 



Dr. 



2An-l + 2X'^An-l + 2}?Bn-l + 2}?An-l + 2\^Bn-l — 2\^Cn-l — 2A^C„_i 

+\''An-l + A^S„_i + X'^Bn-l + X'^Dn-l (H5) 

2 A An~l + 2 A Bn-l + 2 A Bn-l + 2yDn-i + 2An-l + 2Bn-l + 2An-\ + 2Bn-\ 

—2Cn-\ — 2Cn-\ (H6) 

A^A„_i + y^Bn-Y + X^Bn-1 + X^Dn-1 + A„_i + Bn-1 + An-1 + -Bn-1 

— Cn-1 — Cn-1 (H7) 

2A2A„_i + 2X'^Bn-i + 2X^Bn^i + 2X^Dn-i, (H8) 



44 



where we assigned A^'' to the valence-bonds on the n-th cell as the contribution from X{U)\{U'). 
Here N^, is the number of the bonds which share a p-site. The corresponding matrix is 



/2 + 4A2 + A4 4A2 + 2A4 

4 + 2A2 4 + 4A2 

2 + A2 2 + 2A2 
V 2 4 

The initial and final vectors are / = (^o, Bq, Cq, Dq)'^ = (1, 0, 0, 0)-^ and F = {An, Bn, Cn, 
= (1,2,0,1)^, respectively. We find Bn = 2Cn- Therefore, the matrix is reduced to T„ in ( |4.82D 
with the base {An, Cn, Dn)'^ ■ The corresponding initial and final vectors are / = {Ao,Co, Dq)'^ 
= {l,0,Of and F = {An, Cn,DnV = (1,4,1)^. 



-4A2 A^ \ 

-4 2A2 

-2 A2 

2 1 



(H9) 



2. Expectation value of the number operator 



The geometric representation of the expectation value is (|3.29|), where an example of the graph 
V U V is shown in Fig. ^9|(c). We first derive the matrix N^^'. We have 9 kinds of graphs on 
the i-th ceh (Fig. |5^ (al)-(a4), (a6), (b3), (b8), (cl), and (c2)). Let A'f'^-D^^'^ be the quantities 



.-|^y/2-i ^Yie valence bonds which 



respectively. When we attach Af^^-Dl':\ to 



classified by the conditions (i)-(iv) in Appendix 
the i-th cell, we adopt the following rule to assign weight 
share the p-site. When we attach Cj_i to the i-th cell, we assign coefficient —1. For the degenerate 
loop on Ai and Bi, we assign weight 1. For the valence bonds on Cj, we assign coefficient 1/2 to 
cancel weight 2 which is assinged at another end of the loop where we regard it as the end of a 
loop without the number operator. The recursion relations are 



A. 



(p) 



B} 



a 



X'^Ai-i + }?Bi-i + )?Ai-i + X^Bi-i — }?Ci-i — }?Ci-i + A^Aj_i + X'^Bi^ 
+A^B,;_l + A^A-i 



X^Ai_i + A2s,_i + X^Bi^i + A2 A^i 



A^^) = 0. 



(HIO) 
(Hll) 

(H12) 
(H13) 



The corresponding matrix is 



/ 2A2 + A^ 
A2 



V 







2A2 + 2A4 
2A2 
A2 




-2A2 






A^ 
A2 





(H14) 



From the equality Bn = 2C„, it is reduced to N^^^ in ( 4.83|) with the base {An, Cn, Dn 



We derive the matrix N^'^\ We have 16 kinds of graphs on the i-th cell (Fig. ^ (al)-(a3), (a5), 



(bl), (b2), (b4), (b6), (b7), (blO), and (c3)-(c8)). 



Let A^'^^-D^'^^ be the quantities classified by the 



conditions (i)-(iv) in Appendix HI, respectively. From the same rule in the derivation of N, 
recursion relations are 



(p) 



the 



id) 



B, 



id) 



a 



id) 



D 



id) 



-A Bi^i + -A A-i + 2^*"i 2^*^^ 2^^*^^ ~ 2^*^^ ~ 2^*""^ 



X^Bi^i + A^A-i- 



The corresponding matrix is 



/ 1 + A2 2A2 + iA^ -2A2 iA^ \ 

1 2 + A2 -2 A^ 

1 1 _L i\2 1 1 \2 

V 



1 + iA^ -1 iA^ 



1 



I J 



Prom the equahty 5„ = 2C„, it is reduced to ^ in ( [4.83|) 



(H15) 
(H16) 
(H17) 
(H18) 



(H19) 



3. Spin correlation function 



The geometric representation of the expectation value is (3.18|), where an example of the graph 



y U y is shown in Fig. [49Kd). We evaluate the weight of the graph with the spin operators. Since 
there are two valence bonds in a cell, d{i,j) = \i — j + 1\ + le in ( |3.18| ) where le is the number of 



the bonds which are not on the n-th cell {i < n < j). (For example, /g = 6 in Fig. 4£(d)). We 
use the same rule in Appendix p 1| for C„ and that in |H 3| for the assignment of the coefficient |. 



Accordingly, the coefficient (— i in ( 3.18 ) is automatically taken into account. We derive the 
matrix S„. We have four kinds of graphs on the n-th cell (i + 1 < n < j — 1) (Fig. |5l|(a)-(d)). Let 
Sn be the quantity defined by a sum. The sum is taken over V UV on the lattice A„ such that 



the graph consists of the cells shown in Fig. 50 on the A;-th cell {1 < k < i — 1) and those shown in 



Fig. 51(a)-(d) on the /-th cell (1 < / < n). The recursion relation is 



Sn 



n— 1) 



(H20) 



and we have = —2. 

We derive the matrix sf'^^\ We have two kinds of graphs on the i-th cell (Fig. |5l|(e) and (f)). 
The recursion relation is 



S, 



ip) 



^ X X'^Ai^i + i X X'^B.^i + ^ X X'^Bi^i + i X A^A^i. 



(H21) 



The corresponding matrix is (^A^, A^, ^A^, 0) and is reduced to S^'^^^ in ( 4.86| ) with the base 

{Am Cm Dn)'^ ■ 

We derive the matrix S^'''^^ We have two kinds of graph on the j-th. cell (Fig. ^l](g) and (h)). 
The recursion relation is 



S' 



ip) 



1 
2 



-- X 2X'^Si-i - - X 2X^Si-i. 
2 



(H22) 



4fi 



The recursion relations to the j + 1-th cell are 



B 



2\^S'f + 2S'f + 2S'f 



D 



'i+i 



and the corresponding matrix is 



/ 2 + 4A2 + \ 
4 + 2A2 
2 + A2 
2 



(H23) 
(H24) 
(H25) 
(H26) 

(H27) 



A simple calculation leads that the product S'^-^^^S'^-^^ is reduced to T^-^j^S^-^^ with the base 
{An,Cn,DnY and we obtain sj^'^^^ in {WM )- 

We derive the matrix S^'^'^^ We have six kinds of the graphs on the i-th cell (Fig. |5l|(i)-(n)). 
The recursion relation is 

- 



S) 



^ X Ai^i + ^ X Bi^i - ^ X Cj_i - ^ X d-i. 



(H28) 



The corresponding matrix is (|A^, ^A^, — A^, 0) and is reduced to sf'*'^'' in ( 4.86 ). 

We derive the matrix S^'^'^^ We have six kinds of graphs on the j — 1-th cell (Fig. ^(o)-(t)). 
The recursion relation is 



= 

id) _ 



B 



\ X 25.-2 - ^ X 25,_2 



Cf\ = 0. 



1 c 1 e 

- X - - X 



(H29) 
(H30) 

(H31) 
(H32) 



The corresponding matrix is S'^_^ = (0, —2, —1,0). The recursion relations to the j-th cell are 



Af 



2)?Bj—\ — 2}?Cj—i — }?Cj~\\'^Bj—\ (H33) 



{d) 



Cf 



2}?Bj^\ 2Bj^\ — 2Cj—\ — 2Cj—\ 



A Bj-i + Bj-i — Cj-i — Cj-i 



D 



id) 



2A Bj-i, 



(H34) 
(H35) 
(H36) 



47 



and the corresponding matrix is 



2 + 2A2 -4 
l + A^ -2 
\0 2 0/ 



(H37) 



A simple calculation leads that the product S'^*^^ S'^l^ 
{An,Cn,Dnf. We obtain sf'^'^^ in ( p6|) . 



is reduced to S^'^^S,_i with the base 



4. Singlet-pair correlation function 



The geometric representation of the expectation value is (3.11), where an example of the graph 
V U V' is shown in Figs. |49|(e) and (f) for (4.55)-(i) and (iii), respectively. We derive the matrix 
H„. On the n-th cell we have one kind of graph (Fig. ^2|(a)). Let Hn be the quantity defined by a 
sum. The sum is taken over VUV' on the lattice A„ such that the graph consists of the cells shown 
in Fig. |5^ on the k-th. cell {1 < k < i — 1) and that shown in Fig. ^(a) on the l-th. cell {i <l <n). 
It is represented diagrammatically 



and the recursion relation is 



(H38) 



Hn 



2H. 



n-l- 



(H39) 



Therefore, we obtain H„ = 2. 

We derive the matrix H^''-^^ (We consider the case ( 4.55D -(i).) We have one kind of graph on 
the z-th cell (Fig. p^(b)). The recursion relation is 



H. 



(p) 



i-i- 



(H40) 



The corresponding matrix is A^(l,2, 1,0) and is reduced to H^'^^^ in ( [4.86|) with the base 

We derive the matrix H^^^\ (We consider the case (g!55| )-(i).) We have one kind of graph on 
the fe-th cell (|5^(c)). The recursion relation is 



-2X^Hk-i. 



(H41) 



The recursion relations to the k + l-th cell are 



2X^Hf + X'Hf 



B, 



k+l 



c, 



k+l 



D 



k+l 



2\^H'^t + 2H'^^ + 2H'^^ 



\2tti(p) I tti(p) I tti(p) 



2\^H'f. 



(H42) 
(H43) 
(H44) 
(H45) 
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Therefore, we have W'-^l^ = S'gi- 



A simple calculation leads that the product H'^^^-^H'^^ is 



reduced to H^+f Hfc with the base C„, L>„)^. We obtain H^;^^^ in ( [ilSSD 



We derive the matrix ("We consider the case ( 4.55D -(iii).) We have one kind of graph 

on the i-th. cell (Fig. ^(d)). The transfer matrix is obtained by replacing the graph in ( H40| ) by 



Fig. 11(d). We find H'J*^ = Wf and obtain Hf'^"^^ in ( ^ ). Similarly, we find H' 



R,(dp) 



using the graph Fig. |52|(e) . We obtain H^^^^ in ( [4.88| ) 



k+1 



TT/(P) 
^ fc + 1 



5. Correlation function (Cj^cj^ 



The geometric representation of the expectation value is ( |3.23| ), where an example of the graph 
VUV for {(?^^cPj J) is shown in Fig. 49(g). For a line, we have 12 kinds of graphs on the n-th cell 
(f + 1 < n < j — 1) (Fig. ^). We distinguish four cases with respect to two d-sites at left end of the 
lattice A„: (i) the upper site belongs to a line and is an end of the line; (ii) the lower site belongs 
to a line and is an end of the line; (iii) both sites belong to a line and the lower site is an end of 
the line; (iv) both sites belong to a line and the upper site is an end of the line. Let Kn-Nn be the 
quantity defined by a sum. The sum is taken over V UV on the lattice An such that the graph 
consists of loops shown in Fig. ^ on the k-th cell {1 < k < i — 1) and a line shown in Fig. 53 



on 



the /-th cell (i < k < n). The restriction for the sum is that the n-th cell is classified to (i) for 
(ii) for Ln, (iii) for M„, (iv) for Nn- They are represented diagrammatically 



Nn 



(H46) 
(H47) 
(H48) 
(H49) 



We use the same rule as that in ( 4.37] ), when we attach Kn-i-Nn-i to the n-th cell. The recursion 
relations are 



Kn 



Mn 



Nn 



— —Kn-1 — Nn-1 — Ln-l 

The corresponding matrix is 



-Ln-l — Kn~l — }?Ln-l — A^M„_i — }?Kn-l — X^Nn-l 



Kn-l — Ln-l — A^Kn-i — X^Nn-l — A^L„_i — A^M„_i 
Ln-l - Mn-1 - Kn-l - Nn-1 

Mn-1. 



(H50) 
(H51) 
(H52) 
(H53) 



/1 + A2 1 + A2 A2 A2\ 

1 + A2 1 + A2 A2 A2 

1 111 

VI 1 1 1 y 



(H54) 



From the equalities Kn = Ln and M„ = Nn-, it is reduced to G„ in ( [4.901 ). 

We derive Gf'^^\ We have eight kinds of graphs on the i-th cell (Fig. |5^(pkl)-(pm6)). Let 
K^^'' and M'f^ be the quantity defined by a sum. The sum is taken over V VJV on the lattice A, 
such that the graph consists of the cells shown in Fig. on the A:-th cell (1 < k < i — 1). The 



4Q 



restriction for the sum is that the i-th cell is represented by Fig. |5^(pkl) or (pk2) for Kn and by 
(pml)-(pm6) for They are represented diagrammatically 



(p) 
(p) 



(H55) 
(H56) 



We use the same rule described above ( [4.37 ) except for that for Cj. When we attach the graph Ci_i 
to the i-th cell, We assign one more coefficient —1 to cancel the coefficient +1 at the beginning of 
the line where we regarded the right end of the line as that of a loop. The recursion relations are 



(p) 



-AA 



i-l 



M. 



ip) 



= -AAi_i 
The corresponding matrix is 



A-Bj-i — AAj_i — XBi^i + XCi-i + XCi-i. 



A 2A A 
2A 2A -2A 



(H57) 
(H58) 

(H59) 



and is reduced to G^'^^'* in ( 4.9C| ) with the base {An,Cn, DnY^ 



We derive G 



L,{P) 



We have ten kinds of graphs on the j-th cell (Fig. 



^(al)-(c4)). Let Af\ 

Bj^'^ , and cj^^ be the quantity defined by the right-hand side of ( 3.23| ) on the lattice Aj . The sum is 



taken over the graph V UV such that the graph consists of loops shown in Fig. [5^ on the k-th. cell 
(1 < A; < i — 1) and a line shown in Fig. ^ (gl)-(gl2) on the l-th cell {i < I < j). The restriction 
for the sum is that the j-th cell is represented by either Fig. |5^(pal) or (pa2) for A^J'\ (pbl)-(pb4) 



for Bf, and (pcl)-(pc4) for cf. 



They are represented diagrammatically 



A 



B 

C. 



(p) 
j 

(p) 
i 

(p) 



(H60) 
(H61) 
(H62) 



When we attach the j — l-th graph to the j-th cell in C^p\ we assign coefficient 1/2 to cancel weight 
2 which is assinged at end of the line where we regard it as the left end of a loop. The recursion 
relations are 



i(p) 



-XK 



B 



ip) 



-Air,_i 



XK 



AMj_i - XKj^i 



XMj^i 



C 



ip) 



1am,_ 



D 



ip) 



(H63) 
(H64) 

(H65) 
(H66) 



The corresponding matrix is 



f 2X \ 
2A 2A 
A A 
V / 



(H67) 



and is reduced to G^'^^'' in 



We derive . Let A\ -Cl be the quantity defined by a sum. The sum is taken over 

V UV on the lattice Aj such that the graph consists of the cells shown in Fig. ^ on the k-ih. 
cell {'i < k < i — 1). The restriction for the sum is that the i-th cell is classified by the condition 



(i) for A^'f'^ (ii) for B^""' , and (iii) for Cf"' in Appendix HI, respectively. They are represented 
diagrammatically 



A) 



id) 



B 

C. 



(d) 

i 

id) 



After a simple calculation we find A^"^^ = Ai, B, 



Bi, C 



(d) 



(H68) 
(H69) 
(H70) 

Ci, and Di = 0. We use the 



same rule for the derivation of G"' . Using the graph shown in Fig. ^(dkl)-(dn2)), the recursion 
relations to the i + 1-th cell are 



L 



Mi+i 
Ni+i 



-X^Af^ -X-'Bf -Cl"^ -Af 



(H71) 
(H72) 
(H73) 
(H74) 



The corresponding matrix is 



/1 + A2 A2 A2 0\ 
1 + A2 A2 A2 
1 1-10 

VI 1 -1 oy 



(H75) 



From the equalities -ftTj+i = -Li+i and Mj+i = A'j+i, we obtain Gf'^"'^ in (^^ ). 

We derive G^'^^\ Let K^'^\ Lj^\ and ivj'^^ be the quantity defined by the right-hand side of 
( |3.24 ) on the lattice Aj. The sum is taken over the graph V UV such that the graph consists of 
loops shown in Fig. 50 on the k-th cell {1 < k < i — and a line shown in Fig. |5^ (gl)-(gl2) on the 
l-th cell {i < I < j). The restriction for the sum is that the j-th cell is classified by the condition 



(i) for -ftTj-^^ (ii) for L^"'' , and (iv) for N^""' in Appendix |H 5| , respectively. They are represented 
diagrammatically 



-(d) 



Ad) 



N. 



id) 



(H76) 
(H77) 
(H78) 



We use the same rule in the derivation of G^'''^^ Using the graph shown in Fig. |53|(dal)-(ddl)), 



the recursion relations are 



A 



2X^Kf 



X^Lf +X^Kf + X''Nf^ 



(H79) 



.^1 



B 



Cj+i 



D 



The corresponding matrix is 



2Kf^ + 2Nf\ 



( 2A2 + -A2 A^ ^^ 

2 + 2A2 -1 2A2 

1 + A2 -i A2 

V 2 2 / 



(H80) 
(H81) 
(H82) 



(H83) 



From the equahty i?„ = 2C„, the matrix is r educed to G^,'''^^ in (pjOQ with the base C„, A 
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FIGURES 



FIG. 1. The correspondence between the sites x (g An) and r (g C„). 

FIG. 2. Diagrammatic representation of the valence bonds. The sohd (broken) hne represents 
a valence bond in the ket (bra). We distinguish four kinds of bonds and call them as non-closed 
"ket-bond" (a), non-closed "bra-bond" (b), self-closed "kct-bond" (c), and self-closed "bra-bond" 
(d). A d-site with ?7 = oo is denoted by a solid circle, and a p-site with ?7 = is denoted by a 
circle. 

FIG. 3. Examples of valence-bond configurations (a) V, (b) V, and (c) their overlap V U V'. 
The lattice is constructed from a cell with one d-site and four p-sites. P-sites in adjacent cells 
actually is a single p-site. The thin broken lines represent hoppings of electrons. The graph VUV 
is factorized into four connected subgraphs (d)-(g). 

FIG. 4. Diagrammatic representation of the identity h\. yh]y^^ = —bj^ ybl. .,. 

FIG. 5. Examples of the elimination of a p-site with four bonds, which exhaust all the cases. 
For (a), we apply the identity bl.^ybj^^^ = —b]j yb\.^^ once. For (b) and (c), we apply it two times. 

FIG. 6. Example of the procedure of the elimination of the p-sitcs with four bonds. We have 
self-closed bonds (a), degenerate loop (b), and non-degenerate loop (c) after the procedure. 

FIG. 7. Example of graph where sites x and y belong to a single subgraph. 

FIG. 8. (a), (b) Type (i) configuration, (c), (d) Type (ii) configuration. 

FIG. 9. Diagrammatic representation of the identity cl.^^bl. y = ~c^,ct^,j/ (d)- A pentagon 
represents the operator cl^^. 

FIG. 10. (a) and (b) Examples of the procedure of the elimination of type (ii) configuration. 

(c) Example of a line. 

FIG. 11. (a) Type (i) configuration, (b) Type (ii) configuration. 

FIG. 12. (a) Type (i) configuration, (b) Type (ii) configuration, (c) Type (iii) configuration. 

(d) Type (iv) configuration. 

FIG. 13. Model A. (a) A cell compose of two d-sites and one p-site. (b) The lattice constructed 
from the cell with cell labelings. (c) The same lattice as (b) drawn differently with unit cell 
labelings. 



FIG. 14. The dispersion relations £'+ and E_. (a) The parameters are Ai = A2 = 1. (b) The 
parameters are Ai = A2 = 

FIG. 15. Examples of configuration of the valence-bonds in the geometric representation for 
(a) the ground state, (b) the norm of the ground state, (c) the expectation value of the number 
operator n^^, (d) the density correlation function for p-sites, (e) the correlation functions {(^^c^l^) 

and (f) the correlation function {cf^^c^^^). 

FIG. 16. Configurations of the valence-bonds on a cell. The weights including the contributions 
from \{U)X{U') are (a) 2X1, (b) 1, (c) 2A?Ai, and (d) 2Af . 

FIG. 17. Configurations of the valence-bonds on the i-th cell in the geometric representation 
for {n^^^) (a)-(c) and (nf^) (d) and (c). The weights including the contributions from X{U)X{U') 
are (a) Ai, (b) 1, (c) A?, '(d) Ai, (e) AfAi, and (f) A?. 

FIG. 18. Occupation on the p-site (a) and the d-site (b) for A2 = 0.1, 0.5, 1, 2, and 4 (solid 
line) and A2 = Ai (broken line). 

FIG. 19. Density correlation function for the nearest-neighbor p-sites for A2 = 0.5, 1, 2, 4, and 
8 (solid line) and A2 = Ai (broken line). 

FIG. 20. Density correlation function for the nearest-neighbor d-sites for A2 = 0.5, 1, 2, and 4 
(solid line) and A2 = Ai (broken line). 

FIG. 21. Configurations of the valence-bonds on the n-th cell (a) {i + l<n<j — \) in the 
geometric representation for (cf^cj^). Those on the i-th cell (b) and (c), and the j'-th cell (d) and 
(e) for a = p, and on the z-th cell (f) and the j — 1-th cell (g) for a = d. The contributions from 
X{U)X{U') are (a) A1A2, (b) Ai, (c) AiAi, (d) A2, (e) Af A2, (f) A1A2, and (g) A1A2. 

FIG. 22. Momentum distribution functions for the p-site (a) and the d-site for (b) for Ai = A2 = 
0.1, 0.2, 0.5, 1, 2, and 100. 

FIG. 23. Momentum distribution functions for the p-site (a) and the d-site (b) for A2 = 0.01, 
0.1, 0.5, 1, 2, and 100 with Ai = 1. 

FIG. 24. Correlation lengths of the correlation function (Cj^^c^^) (a) and that of the density 
correlation function (b) for Ai = A2. 

FIG. 25. Model B. (a) A cell compose of two p-sites and one ci-site. (b) The lattice constructed 
from the cell with cell labelings. 

FIG. 26. The dispersion relations Ej^ and E^. The parameters are Ai = A2 = 1. 



FIG. 27. Occupation on the p-site (a) and the d-site (b) for A2 = 1, 2, 4, and 8 (sohd hne) and 
A2 = Ai (broken Hne). 

FIG. 28. Density correlation function for the nearest-neighbor p-sites for A2 = 0.5, 1, 2, 4, and 
8 (solid line) and A2 = Ai (broken line). 

FIG. 29. Density correlation function for the nearest-neighbor d-sites for A2 = 0.5, 1, 2, 4, and 
8 (solid line) and A2 = Ai (broken line). 

FIG. 30. Spin correlation functions for the nearest-neighbor p-sites (a) for A2 = 2, 4, and 8 

(sohd line), d-sitcs (b) for A2 = 1, 2, and 4 (solid line), and p- and cZ-sites (c) for A2 = 1, 2, 4, and 
8 (solid line). The broken lines are for A2 = Ai. 

FIG. 31. Singlet-pair correlation function for the nearest-neighbor p-sites for A2 = 1, 2, 4, and 
8 (solid line) and A2 = Ai (broken line). 

FIG. 32. Momentum distribution functions for the p-site (a) and the d-site (b) for Ai = A2 = 
0.01, 1, 10, and 100. 

FIG. 33. Momentum distribution functions for the p-site (a) and the d-site (b) for A2 = 0.01, 
0.1, 1, 5, and 10 with Ai = 1. 

FIG. 34. Correlation lengths of the correlation function (q^c^^) (a), that of the density corre- 
lation function (b) , and that of the spin and the singlet-pair correlation functions (c) for Ai = A2 . 



FIG. 35. Model C. (a) A cell compose of one ^?-sites and four d-site. (b) The lattice constructed 
from the cell with cell labelings. 

FIG. 36. The dispersion relations Ei, E2, and E^. The parameters are (a) A = 1, (b) A = \/2, 
and (c) A = 2. 

FIG. 37. Occupation on the p-site (a) and the d-site (b). 
FIG. 38. Density correlation functions for the nearest-neighbor p-sites (a) and d-sites (b). 

FIG. 39. Spin correlation functions for the nearest-neighbor p-sites (a), p- and d-sites (b), and 
d-sites (c). 

FIG. 40. Singlet-pair correlation function for the nearest-neighbor p-sites. 

FIG. 41. Momentum distribution functions for the p-site (a) and the d-site (b) for A = 0.1, 0.5, 
1, 2, and 10. 



FIG. 42. Correlation lengths of the correlation function (cj^cj^.) (a), that of the density corre- 
lation function (b), and that of the spin and the singlet-pair correlation functions (c). 



FIG. 43. (a) The same lattice as Fig. 25(b) drown differently with unit cell labelings. Examples 
of the configurations of the valence-bonds in the geometric representation for (b) the ground state, 
(c) the norm of the ground state, (d) the norm of the ground state after the elimination of p-sites 
with four bonds, (e) the expectation value of the number operator ra^g., (f) the spin correlation 
function {S^'^Sj'^), (g) the singlet-pair correlation functions {b\jbp,f}i and (h) {h\-hf^j), and (i)the 
correlation function (cf j^c^^^.). (j) Examples of the loops extending over more than two cells. 

FIG. 44. Configurations of the valence-bond on a cell. The contributions from \{U)\{U') are 
(a) Ai, (b) A|, (c) AfAi, (d) Af, and (e) Af. 

FIG. 45. Configurations of the valence-bonds on the z-th cell in the geometric representation 
for (nf^) (a)-(d) and of (nf^) (e)-(g). The contributions from \{U)X{U') are (a) A^, (b) A|, (c) 
A?Ai, (d) Af, (e) Ai, (f) and (g) Af. 

FIG. 46. Configurations of the valence-bonds on the i-th cell (a) and the j — 1-th cell (b) in the 
geometric representation for and on the i-th cell (c) and the j-th cell (d) for {S^^'^Sj''^). 

The contributions from X{U)X{U') are (a) Af A^, (b) Af A^, (c) Af , and (d) A^. 

FIG. 47. Configurations of the valence-bonds on the n-th cell (a) {i + l<n<j — 1) in the 
geometric representation for the singlet-pair correlation functions. Those on the i-th cell (b) and 
the k-th cell (c) and (d) for {bl^bk^k), and on the i-th cell (e) and the k-th cell (f) in {b\jbk,i)- The 
contributions from X{U)\{U') are (a) Af A^, (b) Af A^, (c) Af, (d) Xf, (e) Af A2, and (f) AiA^. 

FIG. 48. Configurations of the valence-bonds on the n-th cell (a)-(c) {i + l<n<j — 1) in the 
geometric representation for (c^g-cjj^). Those on the i-th cell (d)-(f) and the j-th. cell (g)-(j) for 
a = p, and on the i-th cell (k) and (1) and on the j-th cell (m) and (n) for a = d. The contributions 
from X{U)XiU') are (a) A1A2, (b) A?A2, (c) AiAi, (d) A1A2, (e) AiAf, (f) Af A2, (g) Xl (h) AfAi, (i) 
Af, (j) Af, (k) AiAi, (1) Xl (m) Ai, and (n) AfA2. 

FIG. 49. Examples of the configurations of the valence-bonds in the geometric representation 
for (a) the ground state, (b) the norm of the ground state, (c) the expectation value of the number 
operator n^^, (d) the spin correlation function {S^'^ S^^''^) , the singlet-pair correlation functions (e) 

i^i^k^k) and (f) and (g) the correlation function (cf^^c^J,)- 

FIG. 50. Configurations of the valence-bonds on a cell. The contributions from X{U)X{U') are 
1 for the valence bonds which do not belong to a p-site, A^ for the valence bonds which belong to 
a p-site (except for (a6)), and A"^ for (a6). 

FIG. 51. Configurations of the valence-bonds on the n-th cell (a)-(d) (i -|- 1 < n < j — 1) in the 
geometric representation for the spin correlation functions. Those on the i-th cell (e) and (f) and 
on the j-th cell (g) and (h) for {S^'^Sj'^), and on the i-th cell (i)-(n) and on the j-th cell (o)-(t) 

for {S^'''S]'''). 

.^7 



FIG. 52. Configurations of the valence-bonds on the n-th cell (a) (^ + 1 < n < j — 1) in the 
geometric representation for the singlet-pair correlation functions. Those on the i-th cell (b) and 
the k-th cell (c) for {h\^hf^ and on the i-th cell (d) and the u-th cell (e) for {h\ ^h^ j). 

FIG. 53. Configurations of the valence-bonds on the n-th cell (i+1 < n < j — 1) in the geometric 
representation for the correlation function (cj ^^c^^). Those on the z-th cell (pkl)-(pm6) and on the 

j-ih. cell (pal)-(pc4) for and on the i-ih. cell (dkl)-(dn2) and on the j-th cell (dal)-(dd2) 
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TABLES 



TABLE L The cell construction of the models in Refs. ||l],^,^J^,^ . Parameters (a) are those of 
Refs. |]5J2|,^. When we choose Aj's, which are defined in the cell, those in (b), the parameters (c) 
represents those in column (a). The model with boundary condition in the last column satisfy the 
uniqueness condition in Ref. ||5). Here O and P denote open and periodic boundary conditions, 
respectively. The symbol D denotes that the boundary condition depends on the model. 



author 


cell 


dimension 


(a) fb) 


parameters 




(c) 


uniqueness 


U. Brandt 




d>2 


(t) A, : 


= A (j = 1,2,3,4) 




(-A^) 


{d>3) 


and A. Giesekus 




d>2 


T -^J 


= ^ (i = 1,2,3,4), A5 


= V2X 


-2A 


0,P 


A. Mielke 


see text 


d > 1 




= A 




(-A') 


D 


R. Strack 




d = 1 


- Ai 


= A2 = 1, A3 = Xi, A4 = 


-Xi 


A 









d = l 


T Ai 


= A, A2 = —A, A3 = 1 




1 

A 









d = 2 


T Ai 


= A, A2 = —A, A3 = 1 




1 
A 


0,P 


H. Tasaki 


d-sites 


d > 1 


Aj Aj 






A, 


D 


H. Tasaki d- 


and p-sites 


d > 1 


Aj Aj 






A, 


D 


Model A 




d = 1 


Ai, 


A2, A3 = 1 






o,p 


Model B 




d = 1 


Ai, 


A2, A3 = 1 









Model C 




d = 1 


A. 


= 1 (j = 1,2,3,4), A5 = 


X 




0,P 



